UNCUSSIFIED 

SfeduRiTY  CLA{5l^igATI6KrBT  T^^IS  ^A~dr 


REPORT  DOCUMENTATION  PAGE 


II 

ul 

2 


AD-A208  664 


4.  PERFORMING  ORGANIZATION  REPi 


r*~r  f* 

-  riiii.  ..fi  Wpq 


% 


V 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Unlv  of  Maryland 


e 


IFFICE  SYMBOL 
(If  applkabh)  - 


1b.  RESTRlQiVE  MARKINGS-  r  ii' 
ifil  t  tb 


I  »r, 

rtit 


-^®rr 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  Is  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER<S) 

AFOtBR-TK»  »  9  "  U  7  1  0 


7a.  NAME  OF  MONITORING  ORGANIZATION 

AFOSR/NP 


6c  ADDRESS  (Oty,  Sutt,  and  ZIP  Coda) 

College  Park,  MD  20742 


7b.  ADDRESS  (Ot)'.  and  ZIP  Coda) 

Building  410,  Bolling  AFB  DC 
20332-6448 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

arngp _ 


8b.  OFFICE  SYMBOL 
Of  appikabla) 

_il2 _ 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


AFn<;p-ftR-nnAR 


8c  ADDRESS  raty.  State,  and  ZIP  Coda) 


10.  SOURCE  OF  FUNDING  NUMBERS 


Building  410,  Bolling  AFB  DC 

20332-6448 

PROGRAM 
ELEMENT  NO.  ' 

61102F 

PROJEa 

NO. 

2301 

TASK 

NO. 

A8 

WORK  UNIT 
ACCESSION  NO. 

11.  iiiLt  vrKiuaa  iacunty  aastideation)  t>j  r  *'  w., 

(U)  PROPAGATION  OF  SHORT  WmiOH,  HIGH  POWER  MICR0WAVE-4W 
NEUTRAL  AND  IONIZED  MEDIA 

w 

12.  PfRSCMAL  /^UTHOR(S} 

Dr  Destler  (wuHam  W.)  &  Dr  Striffler  (Charles  D.) 

I 

13a.  TYPE  OF  REPORT 
Final 


|l3b.  TIME  COVERED 

1  Dec  85ro  30  Nov 


14, 

L 


DATE  OF 


March 


pS.  PAGE  COUNT 
_ ^ _ 


16.  SUPPLEMENTARY  NOTATION 


17. 


COSATI  COOES 


7 


FIELD 


z:- 


GROUP 


SUB-GROUP 


18.  SUBJEa  TERMS  {Condnua  on  rayana  if  nacassary  and  Idantify  by  block  number) 


4 ')9.  ABSTRACT  (continue  on  ravarsa  If  nacasury  and  idantify  by  block  number) 

>Experiinent8  were  completed  on  the  propagation  of  very  high  power  [(500  kw/cmsq),  short  burst 
(3-30  ms)  microwaves  through  a  neutral  media.  Results  indicate  breakdown  field  strength^ 
are  about  a  factor  of  two  higher  when  the  test  cell  is  shielded,  from  ac^ays  generated  by  the 
microwave  source.  Standard  theoretical  models  were  shown  to  give  very  accurate  predictions 
of  actual  observations  for  all  cases  except  when  the  RF  frequency  was  near  the  plasmJ^'»'  • 
frequency .  Ifeq  W  (A  H.. 

C^cfrOnfc  "tVrCjslioJJd  •  s/»/ Microu)ov'c 

Lac-D^e^.  DrV;V  Ai'r, 


20.  DISTRIBUTION /AVAILABIUTY  OF  ABSTRAO 
XDuNCLASSIFIEDAJNLIMITED  X  □  SAME  AS  RPT. 
22*.  NAME  - - -  -  -  ■ 


1 21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


|22b.  TELEPHONE  f/nc/ude  Area  Coda)  j22c.  OFFICE  SYMBOL 

I - -  1-  ni  wwi.#  i.P 

63  APR  tdttion  may  bt  usad  until  axhaustad. 

Alloth.r.ditions.r.ob»ol«..  - LECURITY  CLASSIFICATION  OF_I  HIS  PAGE 

CLASSIFIED 


8  9  6  0  6  07 §' 


Final  Technical  Report 


AJ'UiMt*  lit*  6  V  *  4^  7  1 0 

STUDIES  OF  THE  PROPAGATIOxN  OF  SHORT  BURST,  HIGH  POWER 
MICROWAVE  RADIATION  THROUGH  NEUTRAL  AND  IONIZED  MEDIA 


Grant  No.  AFOSR-86-0046 


For  the  Period  December  1,  198^ 


to  November  30,  1988 


Submitted  to 

Air  Force  Office  of  Scientific  Research 


Prepared  by 

Electrical  Engineering  Department 
University  of  Maryland,  College  Park,  MD  20742 


UNIVERSITY  OF  MARYLAND 

LABORATORY  FOR  PLASMA  RESEARCH 
COLLEGE  PARK,  MARYLAND 
20742-3511 


Final  Technical  Report 


STUDIES  OF  THE  PROPAGATION  OF  SHORT  BURST,  HIGH  POWER 
MICROWAVE  RADIATION  THROUGH  NEUTRAL  AND  IONIZED  MEDIA 


Grant  No.  AFOSR-86-0046 


For  the  Period  December  1.  1987  to  November  30,  1988 


Submitted  to 

Air  Force  Office  of  Scientific  Research 


Prepared  by 

Electrical  Engineering  Department 
University  of  Maryland,  College  Park,  MD  20742 


Accesion  For  j 

NTIS  CRA&I 

DTtC  TAB 

O 

Unannounced 

0 

1  JliStiflC.ltiOtl 

1  .  1 

By _ _ 

Dist.  ibution/ 


AvdiiaN'ty  Codes 


Di<^t 


Avail  j  'dior 
SiHtctal 


Annual  Tecfanical  Report 


Submitted  to:  Air  Force  OflBce  of  Scientific  Research 

Submitted  by:  Electrical  Engineering  Department 

University  of  Marylauid 
College  Park,  Maurylemd  20742 

Grant  Number:  AFOSR-86-0046 

Period:  December  1,  1987  to  November  30,  1988 

Principal  Investigators:  William  W.  Destler,  Professor 

Charles  D.  Striffler,  Professor 

Title  of  Research  Project:  “Studies  of  the  Propagation  of  Short 

Burst,  High  Power  Microwave  Radiation 
through  Neutral  and  Ionized  Media.” 


2 


I.  INTRODUCTION 


This  Technical  Report  contains  a  summary  of  progress  made  in  AFOSR  funded  studies 
of  short  burst  microwave  propagation  through  neutral  and  ionized  media  over  the  time 
period  December  1,  1987  to  November  30,  1988.  Work  under  this  program  commenced  in 
December  1985. 

The  preliminary  focus  of  this  work  has  been  the  experimental  and  theoretical  study  of 
the  ability  of  very  short  duration  microwave  pulses  to  propagate  through  the  atmosphere 
without  inducing  breakdown.  Significant  progress  has  been  made  in  this  airea  both  in  our 
theoretical  understanding  of  the  breakdown  process  and  in  actual  experimental  studies  of 
air  breakdown  due  to  short  pulse  microwaves. 

A  second  project  supported  under  a  supplementary  greuit  involves  the  study  of  energy 
recovery  in  Large  Orbit  Gyrotrons.  In  this  study,  attempts  are  underway  to  recover  a 
significant  fraction  of  the  energy  of  the  electrons  in  a  Large  Orbit  Gyrotron  zifter  the 
radiation  producing  interaction  is  complete.  In  this  manner,  the  overall  efficiency  of  such 
devices  could  be  dramatically  increased.  Significant  new  experimental  studies  of  such 
systems  have  been  completed  during  the  past  grant  period. 

This  report  is  organized  sis  follows:  A  brief  sununary  of  the  highhghts  of  these  research 
programs  is  contained  in  section  II.  Copies  of  papers  published  during  this  period  sire 
contained  in  the  Appendix. 
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II.  Progress  Report  for  the  Period  December  1,  1987  to  November  30,  1988 

A.  Experimental  Studies  of  the  Propagation  of  High  Power,  Short  Burst 
.  Microwaves  Through  the  Atmosphere 

Initial  experimental  results  of  the  propagation  of  very  high  power  (1-500  kW/cm^), 
short  burst  (3-30  ns)  microwaves  through  a  neutral  media  were  published  in  the  Journal  of 
Applied  Physics  in  a  paper  (“Short-Pulse,  High  Power  Microwave  Propagation  in  the  At¬ 
mosphere”)  enclosed  in  the  Appendix.  During  the  past  year,  we  have  conducted  additional 
experiments  to  determine  whether  x-rays  from  the  microwave  source  could  be  a  source  of 
significant  free  electrons  in  the  test  cell  and  thus  result  in  artificizdly  low  brezikdown  limits. 
These  experiments,  in  which  the  test  cell  was  carefully  shielded  from  the  microwave  source 
by  a  concrete  and  lead  shield  wall,  are  summairized  in  a  second  paper  ( “Shielded  Source, 
Short-Pulse  Microwave  Propagation  Experiments”)  submitted  to  the  Journal  of  Applied 
Physics  for  publication  and  enclosed  in  the  Appendix.  Results  indicate  that  breakdown 
field  strengths  are  about  a  factor  of  two  higher  when  the  test  cell  is  shielded  from  x-rays 
generated  by  the  microwave  source. 

B.  Theoretical  Studies  of  Short  Burst,  High  Power  Microwave  Propagation 

We  have  enclosed  in  the  Appendix  a  copy  of  the  poster  presentation  at  the  1988  APS 
Plasma  Physics  Meeting.  This  work  is  a  summary  of  the  theoretical  models  and  cal¬ 
culations  that  have  been  studied  under  the  existing  contract.  In  our  initial  studies  we 
considered  a  parameter  regime  where  the  ionization  of  the  background  gas  by  the  RF  is 
so  slow  that  the  ionization  rate  and  RF  propagation  time  sczdes  are  very  different.  In  this 
case  we  can  solve  for  the  electric  field  via  the  wave  equation  for  an  inhomogeneous  time- 
independent  plasma  density  profile  and  impose  self-consistency  with  the  ionizing  electric 
field.  We  include  besides  the  ionization  process,  the  attachment  and  elastic  scattering 
processes  for  air.  The  buildup  of  plasma  in  a  gas  cell  is  calculated  along  with  reflection, 
absorption,  and  transmission  properties  of  the  incident  RF  wave.  A  criteria  for  defining 
“breakdown”  is  presented  and  results  are  compared  with  experiments.  The  agreement  is 
found  to  be  satisfactory  for  regimes  where  the  rates  of  ionization,  etc.  are  valid.  These 
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results  have  been  summEirized  in  previous  progress  reports. 

This  past  year  we  have  developed  a  “Medium  Ionization  Rate  Model”  for  RF  propaga¬ 
tion  through  a  plasma.  In  the  model  we  still  assume  the  ionization  rate  is  low  compared 
to  the  RF  frequency  but  do  not  require  the  EM  response  time  of  the  plasma  to  be  small. 
That  is,  finite  EM  transient  time  effects  are  taken  into  account.  For  the  case  when  the 
RF  frequency  is  near  the  plasma  frequency  standau-d  methods  of  calculation  are  no  longer 
valid.  We  have  developed  a  power  series  method  about  the  critical  density  point.  Detaiiled 
convergence  criteria  have  been  checked  and  initial  results  applied  to  experimental  cases 
are  in  progress. 

C.  Studies  of  Energy  Recovery  in  Large  Orbit  Gyrotrons 

These  studies,  fimded  under  a  supplementary  grant  from  AFOSR,  are  summairized  in 
two  papers  enclosed  in  the  Appendix  (“Beam  Conditioning  for  Electron  Energy  Recov¬ 
ery  Systems  in  Devices  Employing  Axis-Encircling  Beams,”  published  in  the  International 
Jovumal  of  Electrorucs,  and  “Study  of  an  Energy  Recovery  System  for  a  Large  Orbit  Gy- 
rotron,”  presented  at  the  1988  International  Electron  Devices  Meeting).  These  studies 
contain  the  first  experimental  demonstration  of  depressed  potential  collection  of  rotating 
beam  electrons  after  their  energy  has  been  reconverted  to  axial  motion.  The  potential  of 
such  recovery  systems  to  recover  70-80%  of  the  energy  remaining  in  the  electron  beam 
after  the  rf  interaction  has  been  clearly  demonstrated  in  these  studies. 
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Beam  conditioning  for  electron  energy  recovery  systems 
in  devices  employing  axis-encircling  beams 

A.  SINGHt.  W.  LAWSONt,  D.  GOUTOSt,  W.  R.  HlXt, 

C.  D.  STRIFFLERt.  V.  L.  GRANATSTEINt 
and  W.  W.  DESTLERt 

Depressed  collectors  are  extensively  used  in  linear-beam  devices.  In  spiralling 
beam  devices,  the  motion  of  the  spent  beam  needs  to  be  converted  into  axial 
motion  before  it  can  deliver  energy  to  the  retarding  axial  electrostatic  held  produc¬ 
ed  by  depressed  collectors.  A  method  of  doing  so  is  presented  here,  along  with 
simulation  results.  A  large-orbit  gyrotron  is  chosen  as  an  illustrative  example  in 
which  the  axis-encircling  beam  is  generated  by  passing  an  axially  streaming  beam 
through  a  magnetic  cusp.  After  rf  interaction,  an  ‘unwinding’  of  the  beam  is  done 
by  a  second  magnetic  cusp  with  a  reversed  magnetic  held.  The  trajectories  of  elec¬ 
trons  going  through  the  two  cusps  have  been  simulated  in  order  to  evaluate  the 
effects  of  hnite  cusp  width  and  hnite  width  of  emitting  surface.  Also  studied  are  the 
effects  of  variations  of  the  magnetic  held,  the  accelerating  voltage  of  the  beam,  and 
the  distance  between  cusps.  These  simulations  show  that  the  rotational  energy  of 
electrons  is  larpMy  converted  into  axial  energy  after  passage  through  the  second 
cusp.  An  initial  geometry  has  been  chosen  for  three  collectors,  of  which  two  are  at 
depressed  potentials.  Changes  in  the  energies  of  the  particles  are  introduced 
between  the  two  cusps  as  an  approximation  of  the  effect  of  rf  interaction.  Simula¬ 
tion  of  trajectories  in  the  region  of  the  depressed  collectors  has  demonstrated 
energy  sorting  of  the  beamlets  to  different  depressed  collectors. 


1.  Introduction 

Gyrotron  oscillators  operating  at  the  fundamental  cyclotron  frequency  have 
made  great  strides  in  achieving  high  rf  power  levels  at  high  frequencies  with  appre¬ 
ciable  efficiencies  of  operation  (Granatstein  1984,  1987,  Fix  et  al.  1984,  Rager  1986). 
Other  device  configurations  using  spiralling  electron  beams  have  unfortunately  suf¬ 
fered  from  significantly  lower  electronic  efficiencies. 

Among  such  devices  are  those  that  use  axis-encircling  beams.  These  include 
large-orbit  gyrotrons  (Lawson  et  al.  1985)  and  circular  geometry  free-electron  lasers 
(Befeki  et  al  1985).  Those  of  the  former  category  of  devices,  in  which  the  axis- 
encircling  beam  is  generated  by  a  magnetic  cusp,  have  been  referred  to  as  cusptrons 
(Uhm  et  al.  1984).  Those  in  which  the  electromagnetic  field  is  supported  by  a 
magnetron-like  structure  have  also  been  called  gyromagnetrons  (Lau  and  Barnett 
1982.  Grow  and  Shrivastava  1982).  The  nature  of  the  beam  and  the  rf  structure 
favours  the  production  of  coherent  radiation  at  a  chosen  high  harmonic  of  the 
cyclotron  frequency.  As  the  value  of  the  magnetic  field  required  for  any  given  fre¬ 
quency  is  reduced  in  the  ratio  of  the  harmonic  number,  this  would  be  a  significant 
advantage  for  operation  at  millimetre  wavelengths  (Chojnack;  al.  1987).  However, 
efficiency  usually  falls  as  the  harmonic  number  rises,  and  this  implies  that  a  con- 
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siderable  proportion  of  the  initial  kinetic  energy  of  the  beam  would  be  carried  by 
the  spent  beam. 

One  possibility  for  improving  the  overall  efficiency  of  such  devices  is  to  recover 
energy  from  the  spent  beam.  In  the  case  of  linear  beam  tubes,  such  as  travelling 
wave  tubes,  this  is  a  well-accepted  approach  (Kosmahl  1982.  True  1987).  The  spent 
beam  is  sorted  into  groups  energy-wise,  and  is  collected  on  electrodes  at  depressed 
potentials.  The  geometry  and  the  potentials  of  the  collectors  are  chosen  to  make  for 
as  soft  a  landing  as  possible  for  each  group.  Thus,  a  good  fraction  of  the  energy  of 
the  beam  is  transferred  back  to  the  power  supply  system.  This  has  the  additional 
advantage  of  reducing  the  bulk  and  cost  of  the  power  supply. 

The  overall  efficiency  of  the  device  tjj  is  related  to  the  intrinsic  efficiency  and 
the  efficiency  of  energy  recovery  by  the  relationship  (Arnush  et  al.  1982) 


hT  = 


>h 

1  -  Ml  -  M 


(1) 


Thus,  if  a  device  of  10%  intrinsic  efficiency  is  operated  in  tandem  with  depressed 
collectors  at  75%  efficiency,  the  overall  elTiciency  would  be  31%. 

In  the  case  of  devices  where  the  spent  beam  has  most  of  its  energy  in  rotational 
motion,  the  electrons  cannot  work  effectively  against  the  electrostatic  field  of  the 
depressed  collectors.  Thus,  it  is  necessary  to  convert  the  rotational  motion  largely 
into  axial  motion  before  leading  the  beam  toward  depressed  collectors.  In  linear- 
beam  tubes,  transverse  motion  in  the  spent  beam  arising  out  of  the  effects  of  the 
focusing  magnetic  field  and  rf  interaction  is  partially  converted  into  axial  motion  by 
an  adiabatic  expansion  of  the  beam  (Kosmahl  1982).  In  the  case  of  axis-encircling 
beams,  the  Larmor  radius  of  the  beam  electrons  is  usually  large,  and  would  be 
further  increased  by  adiabatic  expansion  before  the  process  of  energy  sorting  starts. 
The  latter  pro<.ci5  itself  involves  further  radial  expansion  and  an  increase  of  the 
aperture  of  each  of  successive  depressed  collectors.  Thus,  adiabatic  expansion,  apart 
from  requiring  a  slow  change  of  magnetic  field  over  a  distance,  involves  consider¬ 
able  enlargement  of  the  dimensions  of  axis-encircling  beam  devices. 

We  have  proposed  ‘unwinding’  the  beam  by  passing  it  through  a  second  cusp 
with  a  reversed  magnetic  field.  Its  radial  magnetic  field  component  acts  on  the  azi¬ 
muthal  component  of  velocity  and  converts  it  largely  into  axial  velocity;  this  is  a 
non-adiabatic  process.  As  shown  by  these  studies,  a  substantial  fraction  of  the 
energy  in  the  rotational  motion  of  the  electrons  is  reconverted  into  axial  motion  in  a 
distance  that  is  small  compared  with  the  pitch  of  the  spiralling  motion  and  with 
little  change  in  the  diameter  of  a  hollow  beam. 

The  extent  of  this  reconversion  is  dependent  upon  various  parameters.  We  have 
studied  the  reconversion  phenomenon  as  a  function  of  the  transition  widths  of  two 
cusps,  the  distance  been  them,  the  radial  width  of  the  beam,  the  magnetic  field 
amplitude,  and  the  accelerating  voltage  of  the  beam.  A  single  particle  code  was  used 
to  plot  the  trajectories.  The  results  of  these  simulations  are  presented  here.  They 
verify  the  applicability  of  the  basic  approach  in  that  the  fraction  of  energy  converted 
into  axial  motion  is  close  to  90%  even  with  a  non-optimum  combination  of 
parameters. 

We  have  also  simulated  the  trajectories  of  particles  with  different  energies  in  the 
region  of  the  depressed  collectors.  For  this  purpose  a  system  of  three  collectors,  two 
of  which  are  at  depressed  potentials,  is  chosen  as  an  initial  example.  Energy  sorting 
under  the  combined  action  of  the  magnetic  and  electrostatic  fields  i;;  demonstrated. 
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In  §2,  the  results  of  computer  simulations  relating  to  the  unwinding  of  a  spiral¬ 
ling,  large-orbit  beam  by  a  second  cusp  under  various  conditions  are  reported. 
Section  3  presents  the  trajectories  in  the  collector  region  and  verifies  energy  sorting. 
Section  4  outlines  the  design  considerations  in  the  planned  experiment,  and  conclu¬ 
sions  are  drawn  in  §5. 


2.  Generation  and  unwinding  of  spiralling  electron  beams 

An  axis-encircling  beam  for  a  large  orbit  gyrotron  can  be  generated  by  passing 
an  axially  streaming  beam  through  a  magnetic  cusp,  which  has  a  sharp  reversal  of 
magnetic  field  (Rhee  and  Destler  1974).  The  behaviour  of  an  electron  beam  passing 
through  such  a  cusp  has  been  studied  in  the  cases  of  an  ideal  cusp  as  well  as  a  cusp 
with  a  finite  transition  width  (Scheitrum  and  True  1981,  Rhee  and  Destler  1974). 
The  cases  studied  include  a  balanced  cusp,  which  implies  an  exact  reversal  of  the 
field  on  the  two  sides  of  the  cusp,  as  well  as  more  generalized  changes  of  field  across 
the  cusp  (Destler  and  Rhee  1977). 

A  solution  of  the  equations  of  motion  shows  that  for  a  spiralling  electron 
passing  through  an  ideal  cusp,  the  guiding  centre  radius  upstream  equals  the 
Larmor  radius  downstream  (Rhee  and  Destler  1974).  As  a  special  case  we  note  that 
an  electron  travelling  axially  at  a  radius  r^  through  an  ideal  and  balanced  cusp 
acquires  a  spiralling  motion,  whose  guiding  centre  lies  on  the  axis  and  whose 
Larmor  radius  equals  An  annular  sheet  beam  of  radius  r,,  correspondingly 
emerges  as  a  uniform  rotating  cylinder.  Considering  the  case  of  two  ideal  cusps 
separated  by  a  certain  distance,  the  electron  would  have  a  zero  Larmor  radius 
beyond  the  second  cusp  and  a  guiding  centre  radius  of  Tq,  i.e.,  an  axially-streaming 
beam. 

A  schematic  diagram  of  a  configuration  which  is  the  subject  of  study  here  is 
shown  in  Fig.  1,  A  stream  of  electrons  starts  from  a  ring  cathode  and  is  formed  by  a 
Pierce  gun  into  an  axially-streaming  annular  beam.  After  passing  the  first  cusp,  the 
beam  passes  through  a  region  which  has  an  rf  structure.  The  beam  then  goes 
through  the  second  cusp  (that  has  a  reversed  field),  which  tends  to  unwind  the 
spiralling  motion  before  the  beam  proceeds  toward  the  depressed  collectors. 

In  practice,  the  transition  region  in  the  cusps  has  a  finite  width.  An  approx¬ 
imation  to  the  field  variaii*”’  across  the  cusp  which  has  been  used  with  success  is  as 
follows  (Destler  et  al.  1988.  Lawson  and  Latham  1987): 

B,lro,  z)  =  Bo  lanh  ^  (2) 

c 

where  Bg  is  the  asymptotic  magnitude  of  the  axial  magnetic  field  on  the  two  sides  of 
a  balanced  cusp,  and  (  is  a  parameter  that  is  the  measure  of  the  transition  width 
along  the  z-axis.  In  this  case,  the  guiding  centre  downstream  of  the  cusp  no  longer 
ties  on  the  axis.  Considering  all  electrons  in  a  thin  hollow  beam  starting  at  a  given 
radius,  the  finite  cusp  width  gives  rise  to  a  scalloping  of  the  post  cusp  beam 
envelope.  The  scalloping  wavelength  to  zeroth-order  is  given  by 

2  =  Inroln^  -  1)*'*  (3) 

which  is  the  electron  axial  pitch  distance  for  an  ideal  cusp;  and  where  A  v^^/rgii; 
i),,  is  the  axial  velocity  upstream  of  the  cusp,  D  =  eBo/my  is  the  relativistic  cyclotron 
frequency,  and  y  is  the  relativistic  mass  factor. 
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Figure  I.  Schematic  diagram  of  a  configuration  which  is  the  subject  of  study  here.  A  stream 
of  electrons  starts  from  a  ring  cathode  and  is  formed  by  a  Pierce  gun  into  an  axially- 
streaming  annular  beam.  After  passing  the  first  cusp,  which  causes  the  beam  to  spiral, 
the  beam  passes  through  a  region  which  has  an  rf  structure.  The  beam  then  goes 
through  the  second  cusp  (that  has  a  reversed  field),  which  tends  to  unwind  the  spiral¬ 
ling  motion  before  the  beam  proceeds  towards  the  depressed  collectors. 


For  this  study,  a  parameter  which  acts  as  a  useful  index  of  the  effects  of  the  two 
cusps  is  a.  defined  as  the  ratio  of  the  transverse  to  the  axial  component  of  velocity, 
i.e.. 


I'll 

In  order  to  obtain  some  figures  for  a  typical  set  of  operating  parameters  and  to  gain 
a  better  insight  into  their  interrelationships,  trajectories  were  simulated  under  differ¬ 
ent  conditions.  A  single  particle  code  was  used  for  tracing  the  path  of  the  electrons 
through  the  cusped  magnetic  fields.  The  code  numerically  integrates  particle  orbits 
through  the  cusp  fields.  The  code  uses,  as  a  model  for  the  cusp,  the  axial  component 
of  the  magnetic  field  as  given  in  eqn.  (2).  A  gaussian  profile  is  used  for  the  radial 
component  of  the  field  and  a  Taylor  expansion  of  the  fields  about  r  =  Tq  is  per¬ 
formed.  See  Dcstler  ef  al.  (1988)  and  Lawson  and  Latham  (1987)  for  more  details  on 
the  code.  The  initial  values  of  the  static  parameters  were  close  to  the  experimental 
values  for  the  low  energy  cusptron  operated  earlier  at  the  University  of  Maryland. 
These  were:  Bq  =-  328  8  gauss,  =  1-5 cm,  accelerating  voltage  Fg  =  26 kV.  and 
C  =  4  mm. 

From  such  a  simulation,  an  illustrative  example  of  the  variation  of  a  with  the 
axial  position  is  shown  in  Fig.  2.  The  value  of  a  goes  from  zero  upstream  of  the  cusp 
to  a  finite  value  in  the  region  between  the  cusps,  denoted  by  It  falls  to  a 

low  value  as  the  beam  passes  through  the  second  cusp,  denoted  by  a,„ida,i .  As  the 
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Figure  2.  Illustrative  example  of  the  variation  of  a  =  vjv^  along  the  axial  coordinate.  It 
rises  from  zero  to  a  value  called  as  the  electrons  pass  the  first  cusp,  located  at 

0  cm,  and  falls  again  to  a  low  value  called  after  the  electrons  pass  the  second 

cusp,  located  at  30cm;  =  340  gauss,  r,„„  =  1-5 cm,  Ci. :  =  4mm,  and  Vq  =  26 kV. 

cusps  have  a  finite  width,  the  value  of  turns  out  to  be  finite  rather  than  zero 

as  in  the  ideal  cusp  limit. 

Various  geometric  and  electrical  parameters  influence  the  trajectories,  as 

well  as  These  were  studied  by  varying  the  parameters  one  by  one.  The 

results  are  reported  in  the  following  subsections. 


2.1.  Variation  of  a  with 

For  the  case  of  the  ideal  cusp  where  the  transverse  velocity  component  is  purely 
azimuthal,  the  value  of  is  given  by: 


[cf,  -  {rpSD^r^ 


The  results  of  particle  code  simulations  performed  to  determine  the  effect  of  varia¬ 
tions  of  B„  are  summarized  in  Fig.  3(a).  As  seen  from  eqn.  (5),  a„,j„„,,  is  a  mono- 
tonically  rising  function  of  B^.  It  is  also  seen  that  goes  through  a  cyclic 

variation.  This  variation  is  associated  with  the  radial  scalloping  of  the  electrons 
between  the  two  cusps  and  the  resulting  variation  in  particle  angle  as  it  enters  the 
second  cusp.  The  physical  explanation  of  this  and  cyclic  variations  with  other 
parameters  are  discussed  in  §2.6.  Here  it  may  be  noted  that  as  Bp  increases,  the 
cyclic  variation  occurs  at  shorter  intervals,  due  to  a  progressive  decrease  in  the 
scalloping  wavelength. 

As  a  general  trend,  the  value  of  is  less  than  by  a  factor  of  ten  or 

more.  The  maximum  value  of  is  mostly  below  0  3.  Consequently,  the  energy 

in  rotational  motion  is  less  than  10%  of  the  energy  in  axial  motion. 

For  most  of  the  other  simulation  runs,  the  distance  between  the  two  cusps, 
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Figure  3.  Results  of  particle  code  simulations  relating  to  the  effect  of  variations  in  different 
parameters  on  and  Median  values  of  the  parameters  are:  =  26 kV. 

second  cusp  at  30cm.  r,„„  =  I  5cm.  (,  j  =  4mm.  and  Bo  =  328  Sgauss.  (a)  Effect  of 
variation  of  B„:  (6)  effect  of  variation  of  r,„„,  the  starting  radius  of  the  particle:  (c) 
effect  of  variation  of  the  accelerating  voltage;  and  Id)  effect  of  variation  of  distance 
between  cusps. 
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was  chosen  as  30cm.  For  this  value  of  separation  and  Bo  =  328  8  gauss,  the  value  of 
^'residual  was  0  014.  Clearly,  this  is  the  optimum  position  to  produce  low  transverse 
energy,  and  thus  increased  depressed  collector  efficiency.  This  is  subject  to  a  narrow 
range  of  variation  of  the  starting  radius,  as  discussed  in  the  next  section. 


2.2.  Variation  of  as  with  starting  position  of  electrons 

For  a  given  value  of  Sq,  the  radius  at  which  the  electron  starts,  denoted  by 
influences  i\,  r,, ,  and  a.  For  an  annular  cathode  with  finite  width,  the  value  of 
^midranae  would  thus  depend  upon  the  starting  position  of  the  electron  (substitution 
of  for  Tq  in  eqn.  (5)).  For  this  reason,  it  is  advisable  to  limit  the  percentage  of 
variation  in  r^,^„  to  a  small  value.  Also,  the  choice  of  value  of  should  not  be  too 
close  to  the  threshold  value  of  the  magnetic  field  for  which  the  particles  reflect  (Rhee 
and  Destler  1974),  for  otherwise  the  variations  in  a  due  to  the  variation  in  will 
be  too  large. 

Figure  3(h)  shows  the  variations  of  and  with  the  starting  posi¬ 
tions  of  the  particle.  As  expected.  rises  monotonically  with  r,„„.  Further¬ 
more,  the  value  of  a,„,du,,  varies  cyclically  with  Recall  that  as  changes, 

the  number  of  gyrations  completed  between  the  two  cusps  changes.  Because  of  the 
finite  width  of  the  cusps,  the  beam  would  have  a  scalloping  motion  whose  phase  of 
arrival  at  the  second  cusp  would  depend  upon  resulting  in  the  cyclic  variation 

^feiidual  * 

The  highest  value  of  .  corresponding  to  r,„„  =  I  46cm,  is  0-37.  For  the 

range  of  values  of  for  the  actual  cathode,  fromi  4-1  6cm,  the  lies  for  the 

most  part  between  0  01  and  0-2,  indicating  a  good  percentage  of  energy  conversion 
from  rotational  to  axial  motion.  It  is  observed  that  if  the  distance  between  the  cusps 
is  changed,  the  value  of  for  the  maximum  value  of  ar„,du.i  'S  also  altered. 


2.3.  Variations  of  x's  with  accelerating  voltage 

Trajectory  simulations  were  done  for  a  narrow  range  of  accelerating  voltages. 
Fq,  from  25  5-26-5  kV.  As  seen  in  Fig.  3(c  ),  decreases  slowly  from  2  15  to 

1  94,  while  rises  from  0  (X)5  at  the  median  value  of  V^  to  0-225  and  0  122  at 

the  two  extremes.  These  changes  are  due  to  variation  from  optimum  in  the  phase  of 
scalloping  at  the  second  cusp. 


2.4.  Variation  of  with  distance  between  the  cusps 

As  expected,  the  is  dependent  upon  the  electrical  parameters  and  those 

of  the  first  cusp  only.  Figure  3(</)  shows  the  variation  of  with  the  location  of 

the  second  cusp  with  respect  to  the  first  one.  It  is  seen  that  the  variation  is  cyclic 
between  the  values  0  014  and  0  250.  These  values  are  satisfactory  from  the  point  of 
view  of  energy  conversion.  The  cyclic  nature  of  the  variation  is  discussed  in  section 
2.6. 


2.5.  Variation  of  ot,„,^„,  with  transition  width  of  the  cusps 

As  the  cusp  half-width.  (,  was  varied  from  2  0-5  0  mm  for  both  cusps,  the 
^mtdrifiit  simultaneously  increased  slowly  from  2  02-2  06.  and  the  decreased 
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Figure  4.  Plots  showing  variations  of  with  the  half-width  of  the  second  cusp  going 

from  3-14  mm.  The  parameter  f,,  the  half-width  of  the  first  cusp,  equals  4  mm. 


slowly  from  0  029-0  006.  Thus,  although  an  increase  in  the  cusp  width  increases  the 
shift  of  the  guiding  centre  from  the  axis,  thus  increasing  the  scalloping;  the  changes 
in  a  are  small.  The  values  of  and  depend  essentially  on  Bq-  The  trajectories, 
however,  are  influenced  by  the  modelling  of  the  field  in  the  transition  region.  The 
model  tends  to  be  less  accurate  as  one  moves  away  from  the  median  value  of  r,„„ 
and  near  to  the  corners  of  the  iron  cusp  material. 

Furthermore,  the  value  of  3t„,jdu.i  was  determined  for  r,„„  =  1-46  cm  (the  worst 
case  defined  in  §2.2).  as  well  as  for  other  values  of  r,„„  for  a  variation  of  C  of  the 
second  cusp  from  3-14  mm,  keeping  C  for  the  first  cusp  equal  to  4  mm.  as  shown  in 
Fig.  4. 

It  is  seen  that  for  starting  radii  such  as  l'4Scm,  which  give  a  relatively  high 
®r«>iduji  4t  Cj  =  3  mm,  the  (x„,ida,i  diminishes  as  increases  and  goes  through  a 
broad  minimum  near  Ci  =  10mm.  For  starting  radii  such  as  1-50 cm,  which  gives  a 
very  low  value  of  areiiduxi  for  Cj  =  3  mm,  the  a,„iduii  has  a  sharp  minimum  at  Cj  = 
4  mm  and  progressively  rises  as  Cjincreases  further  to  Nmm. 

The  maximum  swing  in  a,„,dy,i  with  variation  in  r,„„  is  by  a  factor  of  nearly  40. 
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when  both  C,  and  (2  arc  4  mm,  However,  this  factor  is  only  two  when  (2  has  risen  to 
14  mm.  Broadly  speaking,  this  phenomenon  can  be  understood  as  follows.  Changes 
in  occur  due  to  variations  in  r,„„,  because  the  latter  changes  the  phase  of 

scalloping  in  which  the  particle  enters  the  second  cusp.  If  the  second  cusp  is  narrow, 
the  phase  relationship  is  sharply  defined,  resulting  in  clear-cut  reinforcement  or  can¬ 
cellation  of  the  radial  component  of  velocity  due  to  the  two  cusps.  When  the  two 
cusps  are  narrow  and  of  equal  width,  the  cancellation  can  be  nearly  complete,  at  a 
favourable  phase  relationship.  However,  when  the  second  cusp  is  broad,  there  is  a 
correspondingly  large  change  in  the  phase  of  scalloping  during  transit  through  the 
cusp.  Consequently,  reinforcement  or  cancellation  is  not  very  sharply  defined.  The 
general  trend  is  for  to  increase  as  C2  rises. 


2.6.  Cause  of  cyclic  variations  in 

The  solution  of  the  equations  of  motion  shows  that  for  the  case  of  a  finite  tran¬ 
sition  width  of  a  cusp,  the  charged  particle  acquires  a  positive  radial  component  of 
velocity  as  it  passes  through  the  cusp.  This,  in  fact,  causes  a  shift  in  the  positions  of 
the  guiding  centre  away  from  the  axis  and  results  in  scalloping. 

When  the  particle  approaches  the  second  cusp,  its  scalloping  motion  may  have 
given  it  a  positive  or  negative  component  of  radial  velocity,  depending  on  its  phase 
in  the  scalloping  cycle.  If  the  particle  is  moving  outward  as  it  approaches  the  second 
cusp,  the  radial  motion  is  accentuated,  thus  giving  rise  to  a  relatively  larger  value  of 
^rtiKiuii  •  reverse  is  true  for  the  case  of  the  particle  approaching  the  second  cusp 

with  an  inward  motion,  which  then  wholly  or  partly  is  compensated  for  by  the  effect 
of  the  second  cusp. 

In  this  connection,  it  is  relevant  to  look  at  the  trajectories  of  the  particles  from 
an  end  on  position.  Figure  5  shows  such  a  projection  in  the  x-y  plane  for  two 


»(cm)  «(cm) 

(a)  (b) 

Figure  5.  Trajectories  of  the  particles  as  seen  from  an  end-on  position  showing  a  projection 
in  the  transverse  (x-y)  plane  for  two  particles.  In  each  case,  the  larger  circle  represents 
an  end-on  view  of  the  spiralling  motion  in  the  midrange  region  (subscript  1),  and  the 
smaller  circle  represents  the  residual  spiralling  motion  beyond  the  second  cusp 
(subscript  2).  In  (a)  the  particle  is  going  outward  as  it  goes  through  the  second  cusp, 
and  the  circle  showing  residual  motion  is  larger  than  the  case  (b)  where  the  particle  is 
going  inward  as  it  enters  the  second  cusp. 
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Angle  covered  by  iroieclory  uo  to  midpoint  of  second  cusp  (degrees) 

Figure  6.  Plots  of  the  value  of  the  ot„,idy,i  versus  the  transit  an^e,  over  and  above  complete 
cycles  for  the  case  where  the  change  in  the  transit  angle  is  due  to  the  change  in  loca¬ 
tion  of  the  second  cusp.  The  points  represent  snapshots  which  do  not  ail  belong  to  one 
and  the  same  cycle.  1^  =  26  kV,  second  cusp  at  30  cm,  =  I  S  cm,  C,,  j  =  4  mm,  and 
Bo  =  328  8  gauss. 


particles.  In  each  case,  the  larger  circle  represents  an  end-on  view  of  the  spiralling 
motion  in  the  midrange  region  with  a  larger  Larmor  radius.  The  smaller  circle  rep¬ 
resents  the  residual  spiralling  motion  beyond  the  second  cusp,  which  has  a  smaller 
Larmor  radius.  This  second  circle  is  seen  to  have  a  larger  diameter  in  Fig.  S  (a)  when 
the  particle  is  moving  outward  from  the  larger  to  the  smaller  circle  than  in  Fig.  5  (h), 
where  the  trace  connecting  the  two  is  moving  inward.  The  off-axis  locations  of  the 
guiding  centre  are  clearly  seen. 

The  value  of  arctMun  versus  the  transit  angle  (modulo  2n)  was  plotted  for  the  case 
when  the  change  in  the  transit  angle  is  due  to  the  change  in  location  of  the  second 
cusp  and  is  shown  in  Fig.  6.  The  one  cycle  shown  is,  in  fact,  a  synthesis  of  snapshots 
of  points  occurring  over  more  than  one  cycle  seen  in  Fig.  3  (d).  A  similar  curve  was 
obtained  when  the  change  in  transit  angle  was  due  to  a  change  in  the  magnetic  field, 
Bq,  while  keeping  the  distance  between  the  cusps  constant.  This  corroborates  the 
physical  explanation  of  the  cause  of  the  cyclic  variations  discussed  above. 

Correlation  between  transit  angle  and  a„,wu,i  was  also  seen  to  be  good  in  the 
case  of  variations  due  to  a  change  in  C-  However,  the  scatter  of  points  was  larger  in 
the  case  of  variations  caused  by  changes  in  .  This  is  attributed  to  the  fact  that 
when  r,„„  is  changed,  it  changes  not  only  the  and  the  scalloping  wave¬ 

lengths  but  also  shifts  the  trajectories  farther  away  from  the  median  radius  of  the 
hollow  beam.  In  that  case,  the  field  profiles  and  the  effective  cusp  width  vary.  The 
°‘rciiduai  changes  cyclically  with  transit  angle  but  not  on  an  equally  smooth  curve. 


2.7.  yariation  in  the  range  of  scalloping 

The  range  of  scalloping,  as  defined  by  the  difference  between  the  maximum  and 
minimum  radii  of  the  trajectories,  was  studied  in  the  region  between  the  cusps  and 
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Figure  7.  The  range  of  scalloping,  as  defined  by  the  difference  between  the  maximum  and 
minimum  radii  of  the  trajectories  in  the  region  between  the  cusps  and  beyond  the 
second  cusp.  Shown  is  the  case  of  variations  due  to  distance  between  cusps,  Bg  = 
328-8  gauss,  =  1-5 cm,  Vg  =  26kV,and  Cj.  j  =  4mm. 


beyond  the  second  cusp.  Figure  7  shows  the  case  of  variations  in  the  range  scallop¬ 
ing  due  to  variation  in  the  distance  between  the  cusps.  The  cyclic  variations  in 
scalloping  beyond  the  second  cusp  are  seen  to  resemble  the  variations  in  are,Mu.|. 
The  scalloping  in  the  midrange  region  is,  of  course,  dependent  only  on  the  para¬ 
meters  of  the  first  cusp,  and  is  therefore  not  a  function  of  the  distance  between 
cusps. 


2.8.  Correlation  between  analytical  and  numerical  results 

The  parameters  and  can  be  evaluated  both  from  numerical 

simulations  and  the  analytical  relationships  given  in  eqns.  (31  and  (5).  The  results  for 
“midrinic  ^rc  shown  in  Figs.  8  (a)  and  (b)  and  agree  very  closely.  When  the  indepen¬ 
dent  variable  is  Bg ,  the  difference  is  0-4%  for  Bg  =  260  gauss  and  1  %  for  Bg  =  340 
gauss,  with  the  respective  values  of  a  being  10  and  2-5.  Where  the  independent 
variable  is  r,„„ ,  the  difference  is  0-4%  for  the  median  value  of  r,„„  and  rises  to  6% 
for  1  -6  cm.  Plots  of  the  numerical  and  analytical  values  of  A„,„„p  as  a 

function  of  the  starting  radius  are  shown  in  Fig.  9.  The  difference  varies  from  2% 
near  the  median  value  of  r,„„  to  6%  at  r,„„  =  1  -6  cm.  These  differences  are  due  to  a 
combination  of  finite  mean  cusp  width  and  variation  of  cusp  width  about  the  mean 
value. 
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Figure  8.  A  comparison  of  values  of  obtained  from  analytical  and  numerical  tech¬ 

niques.  (a)  Variation  of  <x  with  Bq  ■  The  difference  is  0-4%  at  Bq  »  260  gauss  and  1%  at 
Bo  =  340  gauss;  1^  —  26  kV,  r,„„  =  I -5  cm,  and  Ci  ”  4  mm.  (6)  Variation  of  a  with 
The  difference  is  0  4%  at  =  1  52cm  and  6%  at  1  40  and  1  60cm:  Bo  »■ 
328  8  gauss,  C|. 2  “  4mm, and  Fg  >■  26kV. 
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Figure  9.  A  comparison  of  the  values  of  scalloping  wavelength  versus  starting  radius  of  the 
particle  as  obtained  from  analytical  and  numerical  techniques.  The  difference  varies 
between  2%  and  6%.  Sq  =  328-8  gauss,  Ki  =  26kV,and  =  4  mm. 


3.  Energy  sorting  in  the  collector  region 

With  its  rotational  motion  converted  sufficiently  into  axial  motion,  the  beam 
would  be  ready  to  enter  the  region  of  depressed  collectors.  The  geometry  and  the 
potentials  in  the  collectors,  as  well  as  the  magnetic  field  profile,  have  to  be  suitably 
chosen.  The  objective  is  that  the  combined  effect  of  this  electrostatic  field  and  the 
magnetic  field  cause  the  trajectories  of  electrons  of  different  energy  ranges  to  be 
sorted  out  to  different  electrodes  with  depressed  potentials  as  close  to  the  energy  of 
that  group  as  possible.  The  larger  the  energy  spread  due  to  rf  interaction  and  the 
more  ambitious  the  target  value  for  collector  efficiency,  the  larger  the  number  of 
depressed  collectors  that  would  be  required. 

In  the  present  case,  considering  the  fact  that  a  large-orbit  gyrotron  operating  in 
a  high  harmonic  mode  has  a  low  efficiency,  a  configuration  was  chosen  consisting  of 
three  collectors,  two  of  which  are  at  depressed  potentials.  In  choosing  the  geometry, 
a  consideration  kept  in  mind  is  that  a  small  component  of  radial  electric  field  is 
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introduced  in  the  collector  region  which  opposes  the  motion  of  any  secondary  elec¬ 
trons  emitted  at  the  collector  surfaces  back  toward  the  axis. 

The  Herrmannsfeldt  code  (Herrmannsfeldt  1979)  has  been  used  for  calculating 
the  trajectories  of  the  particles.  The  magnetic  field  profile  is  also  evaluated  by 
specifying  the  position,  the  geometry,  and  the  current  through  each  ideal  N-turn 
coil.  The  effect  of  the  iron  plate  that  generates  the  second  cusp  was  simulated  by  an 
image  of  the  coil  nearest  to  the  iron  plate. 

The  potentials  on  the  depressed  collectors,  as  well  as  the  parameters  of  the  coils, 
were  adjusted  for  suitable  energy  sorting.  The  electrode  potentials  were  set  at  0, 
— 12-5,  and  —  20kV.  The  magnetic  field  was  always  adjusted  to  make  it  a  balanced 
cusp. 

In  fixing  the  parameters  for  trajectory  studies,  the  energy  of  the  electron  beam 
entering  the  first  cusp  was  initially  assumed  to  be  26keV,  corresponding  to  an 
earlier  experiment  on  cusptron  operation  at  the  University  of  Maryland  (Chojnacki 
et  al.  1987).  In  order  to  simulate  the  conversion  of  a  part  of  the  energy  in  the  beam 
into  rf  in  the  region  between  the  cusps,  the  particles  were  divided  into  five  groups. 
Energy  changes  were  introduced  to  correspond  to  -30%,  -20%,  -  10%,  0%,  and 
-1-10%  of  26keV.  The  radial  location  at  which  the  particles  emerged  from  the 
second  cusp,  as  well  as  their  angles,  were  noted  from  single  particle  trajectory  code 
simulation.  These  were  then  used  as  inputs  to  the  Herrmannsfeldt  code  simulation 
of  trajectories  in  the  collector  region.  The  trajectories  arc  shown  in  Fig.  10.  The 
beam  current  was  assumed  to  be  one  ampere,  as  in  the  earlier  experiment. 


f  :(mm)  — 

Location  of 
second  cusp 

Figure  10.  Electron  trajectories  in  the  collector  region  for  electrodes  al  ground.  -  12  5kV, 
and  — 20kV.  The  beamlets  ol  energies  18  2  and  20  8keV  are  collected  at  -12  5kV; 
those  of  energies  23  4.  26  0.  and  28  6kcV  are  collected  at  -20kV. 
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It  is  seen  that  the  beamiets  of  energies  18-2  and  20-8  keV  are  collected  at  the 
electrode  at  —  12  5kV,  and  those  of  energies  23-4.  26  0  and  28-6  keV  are  collected  at 
the  electrode  at  —  20kV.  These  results  show  energy  sorting  in  the  collector  region 
with  two  depressed  collectors  for  a  case  of  small  spread  in  energy,  arising  out  of 
relatively  low  electronic  efficiency. 

For  larger  electronic  efficiency,  it  would  be  necessary  to  include  a  more  precise 
model  of  the  energy  spread  in  the  beam  arising  out  of  the  rf  interaction.  A  larger 
numtier  of  depressed  collectors  would  also  be  desirable  for  the  case  where  electronic 
efficiency  is  high. 

Furthermore,  trajectory  simulation  for  14  5  keV  shows  a  tendency  for  some  elec¬ 
trons  to  be  reflected  before  touching  the  -12  5kV  electrode.  This  indicates  a 
residual  rotational  energy  of  2keV  for  some  electrons.  This  phenomenon  is  mini¬ 
mized  by  suitable  profiling  of  the  magnetic  and  electrostatic  fields.  It  also  points  up 
the  benefit  accruing  from  the  transformation  from  rotational  to  axial  motion,  failing 
which,  for  example,  the  collection  of  the  23  4keV  electrons  on  the  -20kV  electrode 
would  not  have  been  possible. 

Since  the  value  of  the  axial  magnetic  field  is  reduced  as  the  beam  proceeds 
farther  into  the  collector  region,  a  further  conversion  of  the  rotational  motion  into 
axial  motion  occurs.  The  reduction  in  magnetic  field  here  need  not  be  slow  enough 
to  qualify  for  being  called  adiabatic. 

For  comparison,  it  may  be  noted  that  if  the  reduction  in  a  were  to  be  obtained 
purely  by  adiabatic  change  of  magnetic  field,  then 

a  =  [Uo' +  I)/.- 1]'*'^  (6) 

where  i,,  and  s  represent  the  ratio  of  transverse  to  axial  velocity  before  and  after  the 
adiabatic  change,  and  /„  is  the  ratio  of  initial  to  final  value  of  the  magnetic  field. 
Correspondingly,  the  Larmor  radius  of  the  orbit  increases  by  a  factor  of  v/m- 
instance,  starting  with  =  2,  an  adiabatic  reduction  of  magnetic  field  by  a  factor  of 
nine  would  reduce  x  to  0  312  and  increase  the  Larmor  radius  of  the  orbit  by  a  factor 
of  three. 

These  considerations  can  be  used  in  designing  any  particular  system,  depending 
upon  the  nature  of  the  device  and  the  requirements  as  well  as  the  space  that  could 
be  made  available  in  the  radial  and  axial  directions  in  the  collector  region. 


4.  Design  considerations  for  the  experiment 

The  computer  simulation  results  described  in  the  previous  sections  arc  used  to 
provide  guidelines  for  the  design  of  the  experiment.  These  guidelines  are  listed 
below.  Recall  that  there  are  limitations  in  the  modelling  of  the  axial  magnetic  field 
at  the  beam  radius  by  the  function  tanh  z/C,  especially  as  the  particles  travel  to  radii 
that  are  'vdl-removed  from  the  median  radius  of  the  beam. 

I.  As  2  is  a  function  of  starting  radius,  the  range  of  variation  of  the  latter  is  to  be 
restricted  (e  g.  to  5%).  Too  large  a  spread  in  would  reduce  the  effi¬ 

ciency  of  rf  interaction.  Furthermore,  if  r,„„  has  a  narrow  range,  then  the 
other  parameters,  such  as  Bq.  Fq,  and  the  location  of  the  second  cusp,  can  be 
fine  tuned  so  that  the  !x„„dB,i  is  low  for  all  the  trajectories.  However,  as  the 
range  of  r,„„  becomes  narrower,  the  cathode  current  density  for  a  given  beam 
current  would  become  higher;  and  thus  become  a  limiting  factor. 
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2.  In  order  to  minimize  scalloping  after  passage  through  the  first  cusp,  its  tran¬ 
sition  width  is  to  be  minimized.  The  limiting  factor  would  be  the  saturation 
flux  of  the  material  used.  Presently  is  set  at  4  mm,  keeping  the  existing  cusp 
number  one  as  a  starting  point  (Chojnacki  et  al.  1987). 

3.  The  median  values  of  other  parameters  are  chosen  so  that  the  beam  is  scal¬ 
loping  inward  as  it  enters  the  second  cusp,  minimizing  . 

4.  The  choice  of  the  width  of  the  second  cusp  is  a  design  compromise.  When  the 
starting  radius  of  the  particles  is  restricted  to  a  range  less  than  5%  and  the 
other  parameters  finely  tuned,  the  minimum  residual  energy  in  the  beam  as  a 
whole,  associated  with  transverse  motion  beyond  the  second  cusp,  occurs 
when  the  two  cusps  are  of  comparable  width.  Increase  of  the  width  of  the 
second  cusp  increases  the  residual  energy  by  a  few  percentage  points. 
However,  this  gives  the  advantage  of  an  increased  radial  gap  available  for  the 
beam  to  pass  through.  In  this  case,  a  cusp  half  width  of  6  mm  has  been 
chosen. 

5.  The  magnetic  field  profile,  the  electrode  shapes,  and  their  potentials  are 
tailored  to  give  optimum  energy  sorting  between  electrodes  and  soft  landing 
of  electrons  with  minimum  residual  energy. 

6.  The  range  of  collection  of  electrons  from  18  2-21  keV  on  the  -  12-5  kV  elec¬ 
trode  and  22-28  keV  on  the  -20kV  electrode  is  adequate  for  a  device  of 
electronic  efficiency  3:10-15%.  Such  a  typically  low  electronic  efficiency  for 
harmonic  operation  does  not  cause  a  large  spread  in  the  energies. 

We  are  working  toward  an  experimental  confirmation  of  these  ideas.  We  are 
also  considering  geometries  involving  larger  numbers  of  depressed  collectors. 


5.  Conclusions 

The  overall  efficiency  of  devices  employing  spiralling  electron  beams  can  be 
raised  appreciably  by  energy  recovery  from  the  spent  beam.  This  can  be  done  by 
’unwinding’  the  beam  before  leading  it  to  depressed  collectors.  A  scheme  for  doing 
this  has  been  proposed  and  analysed  for  the  case  of  a  large-orbit  gyrotron.  It 
involves  the  use  of  a  second  cusp  with  a  reversed  magnetic  field  for  the  purpose  of 
beam-conditioning.  It  is  seen  that  in  this  way  a  large  fraction  of  the  energy  in  the 
spent  beam  gets  converted  into  energy  of  axial  flow.  Furthermore,  the  following 
conclusions  emerge: 

2rej.<iu.i  varies  slowly  with  accelerating  voltage. 

2-  varies  cyclically  with  starting  radius,  magnetic  field,  and  distance 

between  cusps. 

3.  X, is  a  minimum  when  the  scalloping  motion  of  the  electrons  is  inwards 
as  they  enter  the  second  cusp. 

4.  The  second  cusp  can  be  made  wider  than  the  first  cusp  to  allow  for  a  larger 
gap  to  let  the  spent  electron  beam  through. 

5.  The  combined  action  of  the  electrostatic  field  provided  by  depressed  collec¬ 
tors  of  suitable  shape  and  the  magnetic  field  profile  beyond  the  second  cusp 
sorts  the  electron  stream  coming  out  of  the  latter  onto  the  depressed  collec¬ 
tors. 
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6.  Collection  at  depressed  potentials  results  in  energy  recovery.  This  enhance¬ 
ment  of  overall  efficiency  is  feasible  as  well  as  particularly  attractive  for 
modes  of  operation  that  have  low  electronic  efficiency. 

The  basic  approach  can  be  applied  to  other  types  of  devices  employing  spiralling 
electron  beams,  for  enhancing  their  overall  efficiency. 
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ABSTRACT 

'.:i  approach  towards  energy  recovery  from  the  spent 
til  lit'  a  large-orbit  gyrotron  is  outlined  in  the  form 
-  .1  'louble  magnetic  cusp  arrangement.  The  first  cusp 
"iiirc.-s  a  large  rotational  velocity,  the  second  cusp 
; die  purpose  of  unwinding  the  axis  encircling  beam 
it  proceeds  towards  depressed  collectors.  The  re- 
<  •>!  numerical  simtilations  are  summarized  for  beam 
i.Miiission  through  the  two  cusps. 

\  i-xperimental  set-up  for  proof  of  principle  is  de- 
Results  are  presented  in  the  form  of  depressed 
'■'iitiiil  on  a  Faraday  cup  versus  the  current  collected 
Current  collection  is  observed  up  to  larger  de- 
potentials  when  the  second  cusp  is  introduced, 
■Miiiinrcd  with  the  case  where  field  reversal  is  not  in- 
thus  verifying  the  concept.  A  design  for  mul- 
lii'presiied  collectors  is  presented  along  vdth  trajec- 
■ '  Id  the  collector  region. 

1.  INTRODUCTION 

•  '■  oitjcctive  of  this  work  is  to  demonstrate  energy  te- 
from  the  spent  beam  by  the  use  of  depressed  coUec- 
■  *  di  >piraling  beam  devices  as  a  means  of  enhancing  over- 
*  ■  '’^.'-ii-iicy.  Depressed  collectors  for  linear  beam  deidoes 
‘Mi'n  the  subject  of  extensive  study  and  ^>pfication. 

i 

■  ■<! 'plying  depressed  collector  techiuques  to  spiraliag 
’‘•■'Vires,  it  is  necessary  in  the  latter  case  to  convert 
'">it.mtial  rotational  motion  in  the  spent  beam  into 
‘  ■  '•  !:..ition  before  it  proceeds  to  depressed  collectors.  We 
'  :  :«|>osed  a  scheme  for  studying  this  for  a  large  or- 
‘'•‘‘"'ron,  using  s  double  cusp  arrangement.(5]  Figure  1 
•  ^  M'heinatic  diagram  of  such  a  system.  A  magneUc 

'''  a  reversed  field  profile  Ibllowa  the  rf  interaction 
"'•■n  Aiiti  reconverts  the  rotational  electron  motion  into 
“  **  «i<»iion. 

lin^-e  examined  the  trajectories  of  electrons  going 
'-'ll  die  two  cusp  system.  [G]  In  this  way  the  residual 
motion  left  in  the  electrons  due  to  parameters, 
''  linite  transition  widths  of  the  cusps,  the  distance 
the  magnetic  field  profile,  the  beam  volt- 
'  *'='1  die  radius  at  which  the  electrons  start,  has  been 


estimated.  The  choice  of  parameters  is  in  a  range  defined 
by  the  earlier  ctisptron  experiments  at  the  University  of 
Maryland.  [7]  Illustrative  data  of  these  studies  and  guide¬ 
lines  worked  out  for  optimization  of  the  conversion  process 
are  presented  in  section  2. 

We  have  also  set  up  a  proof-of-prindple  experiment.  Its 
features  are  described  in  section  3.  The  experimental 
results  on  beam  transmission  through  the  two  cusps  and 
measurements  on  the  current  collected  by  a  Faraday  cup  at 
depressed  potentials  are  presented  in  section  4.  It  is  seen 
tb  st  with  the  field  reversal  at  the  second  ertsp,  the  current 
collection  continues  up  to  higher  depressed  potentials  than 
without  such  reversaL 

FVirthermore,  the  electron  trajectories  of  the  spent  beam 
have  been  simulated  in  the  depressed  collector  region,  un¬ 
der  a  range  of  conditions  of  collector  geometry,  potentials, 
magnetic  field  profiles  and  energy  variations  in  the  incoming 
beamkts.  Typical  trajectories  and  estimates  of  collector  ef¬ 
ficiency  are  given  in  section  5,  along  with  the  latest  design. 
Conclusions  are  summarized  at  the  end. 

a.  CONVERSION  or  ROTATIONAL  MOTION  INTO 
AXIAL  MOTION 

b  this  study,  a  useful  parameters  is  a,  defined  as  the 
ratio  of  the  perpendicular  component  of  velocity  Vx  fo 
parallel  component  of  velocity  V|.  A  representative  example 
of  studies  of  OmUmatt  defined  as  vx/o|  in  the  region  after 
the  second  cusp,  is  shown  in  Figure  2,  where  it  is  plotted 
versus  the  starting  radius  of  the  dectrons.  The  starting  ra¬ 
dius  varies  from  1.4  to  1.6  cm,  the  beam  voltage  is  26  kV, 
the  magnitude  of  the  magiaetic  field  b  330  Gauss.  As  the 
particles  cross  the  first  cusp,  the  o  is  a  function  of  the  start¬ 
ing  radius  r„.  The  finite  transition  vridth  of  the  first  cusp 
results  m  the  guiding  center  bring  somewhat  cUsplaced  from 
the  axis.  As  a  result,  the  particle  enters  the  second  cusp 
with  a  radial  motion  that  is  inwards  or  outwards  depend¬ 
ing  upon  the  number  of  gyrations  completed,  especially  the 
fractional  part.  This  fractional  part  converted  into  degrees 
(full  cycle  s  360*),  is  also  shown  in  Figure  2. 

It  is  seen. that  Qrttidaai  has  a  cyclical  variation  with  r^. 
However,  it  stays  generally  below  0.32,  representing  only 
about  10%  of  the  electron  energy  in  residual  rotational  mo¬ 
tion.  The  cyclic  variation  is  linked  with  the  fractional  part 


'  '':H.H,MiWiri(i(Wi|55  $1.00  0  1988  IEEE 


Best  Available  Copy 


lEDM  88-153 


of  the  gyrations  performed  in  transit  between  the  two  cusps. 
As  discussed  in  detail  elsewhere  [6|,  this  can  be  understood 
by  recalling  that  the  finite  width  of  the  first  cusp  results  in 
an  outward  motion  of  an  axially  streaming  particle  passing 
through  it,  so  that  over  and  above  the  conversion  of  the  ax¬ 
ial  to  rotational  velocity,  the  guiding  center  is  shifted 
ofF-axis.  resulting  in  a  finite  OretUumi-  But  a  particle  that 
enters  the  second  cusp  with  an  inward  radial  motion  has  it 
counteracted  by  the  latter,  and  vice  versa.  To  minimize  the 
Ortnduai  over  the  whole  beam,  it  is  desirable  to  reduce  the 
range  of  r„,  consistent  with  the  current  density  available 
from  the  cathode.  The  set  of  parameters  can  then  be  fine 
tuned  to  give  a  low  average  of  The  width  of  the 

first  cusp  needs  to  be  minimized.  That  of  the  second  cusp 
needs  to  be  of  the  same  order,  though  a  somewhat  greater 
value  can  be  accommodated,  so  as  to  provide  a  larger  aper¬ 
ture  for  the  spent  beam  to  pass  through. 

.  3.  THE  EXPERIMENTAL  SET-UP 

The  basic  arrangement  is  shown  schematically  in  Figure 
1.  A  photograph  of  the  set-up  is  shown  in  Figure  3.  In  this 
the  electron-gun  is  on  the  left  hand  side,  followed  by  two 
iron  plates  which  generate  the  magnetic  cusps  with  the  aid 
of  coils  that  carry  appropriate  currents.  There  are  seven  of 
them  in  ail. 

Profiles  of  the  magnetic  field  have  been  generated  to 
correspond  to  near  constant  B,  field  in  the  gun  and  if  re¬ 
gions  denoted  by  Bq,  and  controlled  reductions  of  the  x 
component  in  the  collector  region.  Profiles  have  been  also 
measured  for  individual  contributions  of  each  cotL  If  the 
flux  density  is  well  below  saturation  in  the  iron  plate,  linear 
superposition  can  be  assumed,  and  different  profiles  synthe¬ 
sized. 

In  the  initial  proof  of  principle,  a  Faraday  cup  is  located 
in  the  collector  region.  It  consuls  of  copper  arifh  a  coating 
of  carbon.  This  coating  would  help  to  reduce  the  number 
of  back-scattered  electrons.  It  can  be  put  at  depressed  po¬ 
tentials  by  a  separate  power  supply.  The  pulse  as  well  as  dc 
components  of  the  current  collected  at  different  potentials 
are  measured.  The  wave  form  of  the  current  »  converted 
into  a  voltage  pulse  by  a  pick-up  coil,  wUch  b  then  dis¬ 
played  on  an  oscilloscope.  The  dc  oompooent  of  current  b 
measured  with  the  help  of  a  filter  that  bypasses  the  higher 
frequency  components  of  the  current.  The  cathode  in  the 
gun  is  fed  with  negative  pulses  up  to  30  kV  and  duration 
of  5  microseconds,  at  a  rate  of  60  pulses  per  second.  The 
duty  ratio  b  then  0.0003.  For  a  pulse  of  collector  current 
of  0.2  ampere,  the  dc  component  of  the  collector  current 
would  thus  be  60  microamperes. 

4.  COLLECTOR  CHARACTERISTICS 

Measured  characteristics  of  the  Faraday  cup  that  can 
in  a  limited  way  function  also  as  a  depressed  collector  are 
shown  in  Figure  4  for  one  set  cf  parameters,  corresponding 
to  an  accelerating  voltage  of  15.5  kV  and  a  field  strength 
i7o  of  255  Gauss.  The  cusp  half-width  b  approximately  4 
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mm  for  the  first  cusp  and  8  mm  for  the  second  cusp, 
distance  between  them  is  32  cm.  ^ 

One  curve  is  for  the  case  where  a  magnetic  field  rente 
is  produced  at  the  location  of  the  second  iron  plate,  n 
second  curve  represents  the  case  where  the  direction  of  el 
rent  in  the  coils  in  the  collector  region  is  reversed,  aadfi 
values  so  adjusted  as  to  give  the  same  value  of  Bo  I 
the  second  cusp,  without  reversal  of  direction.  This  u 
an  approximation  to  a  constant  field,  to  contrast  witk  ( 
case  of  field  reversal.  Actually,  there  is  a  dip  in  the  nk 
the  B,  component  at  the  location  of  the  second  iron  phi 
because  the  lines  of  field  tend  to  fan  out  towards  the  ini 
Thb  tends  to  influence  the  beam  transmission  through  k 
hole  in  the  second  iron  plate.  Furthermore,  while  the  la 
face  of  the  cup  b  treated  to  reduce  secondary  emissiaa,| 
does  not  have  a  retarding  field  to  suppress  secondaries,  I 
in  the  case  of  the  design  in  the  next  section.  Consequeath 
the  collector  current  here  is  diminished  by  some  Ktunk 
secondaries.  However,  the  maximum  depressed  potca 
at  which  the  collected  current  reduces  to  zero  still  renaiJ 
valid.  wi* 

It  b  seen  that  the  maximum  depressed  potential  I 
which  the  Faraday  cup  collects  current  b  -13  kV  in  theei^ 
where  the  second  cusp  b  present;  and  -7  kV  to  the  OS 
where  there  b  no  second  cusp.  These  figures  demooiOM 
that  the  introduction  of  the  second  cusp  has  increased  ill 
depressed  potential  up  to  which  the  ‘spent’  beam  caa  I 
collected,  thus  increasing  ertergy  recovery.  The  value  cf  t 
after  a  particle  passes  the  first  cusp  b  pven  by:  , 

where  b  the  axial  vdocity  prior  to  reaching  the  fa 
cusp,  and  w,  b  the  angular  cydotroo  frequency.  The  •  H 
an  accelerating  vdtage  of  15.5  kV,  Bo  —  255  Gauss,  a* 
varying  frocn  1.4  cm  to  1.6  cm,  varies  from  1.57  to 
Thb  implies  that  there  b  a  spread  in  the  energy  of  a^ 
motion,  both  in  the ‘rT  section  and  beyond  the  second  OM 
The  maximum  depressed  potential  of  -13  kV  correspocM^ 
dectrons  which  emerged  from  the  second  cusp  writh  aa  an 
0.43,  indicating  conversion  of  rotational  into  axial  mot^ 
As  discussed  already.  Cor  getting  conditions  closer  to 
the  cusp  widths  as  well  as  the  range  of  would 
be  reduced. 

5.  A  DESIGN  FOR  MULTIPLE  DEPRESSED  ■ 
COLLECTORS 

When  rf  interaction  takes  place,  the  energy  of  the  I 
gets  spread  out  over  a  narrow  or  wide  range, 
upon  the  dectronic  eflSdency.  In  order  to  have  i 

energy  recovery,  it  b  necessary  to  sort  the  beam  into  bea^ 
lets  covering  smaller  ranges  of  energy,  collected  on  multir^ 
cdlecton  each  depressed  to  a  potential  adjusted  Cor  a* 
lecting  the  beamlet  with  minimum  loss  of  energy,  i.e., 
soft-landing.  ^ 

We  have  nm  simulations  of  trajectories  using  the 
mannsfddt  code.  [8j  Systems  employing  three  and 


t56-lEDM  RH 


Best  Available  Copy 


.  •  ••nMthrawl  for  pnitiaiy  ritetron  (riii«cioriN. 
a  t  II  liii  Kiudy.  Aw  raliteton  hftw  bwn  UMit  nnd 
,  <<iii  loUn  (owMdt  cenddmtion  of  the  effect  of 
. ,  11^1 1'lectronf.  Tbeee  aetondery  or  reAected  pri> 

^  >  null  ilieroUectoreAkiency  when  they  end  upon 

,  8  ,ii  .1  \inm  depceeeed  potential.  Tliey  can  diiturb 

trtrimnt  if  tliey  head  all  the  way  ba^  to  that  re- 
M  iiinni  for  the  primary  particlea  are  ahown  in 
.  .UMiiiitfta  are  at  enerxiee  of  18.2,  20.8,  23.4,  20.0. 
.  •  •  kV*.  !he  Matter  repreaenting  the  effect  of  rf  in- 

^  '  8.  iliHYty  forting  ia  obaerved  among  the  deprcaaed 

••  I  • 

•.•^iiK  the  geometry  of  the  collectora  and  valuea  of 
;..>ti«‘im.  a  conatderatlon  kept  in  mind  waa  that  any 
..  i.iiv  I'iiYtrona  or  reAected  primariea  releaaed  from  its 
■  • :  iiiiiMct  should  tee  a  retarding  eleetroitatic  Add.  ao 
...itiiiiue  ion  of  collector  eRideney,  or  backtracking 

•  In  the  caae  of  energetic  reAected  primariea 

•  It  •Miirr  a  region  of  jtceeleration  by  overcoming  the 
'••I  t•■•l,tciing  Aeld.  Thua  the  conAguration  and  param- 

•  -  xitt  fnr  Fig.  5  are  not  thoae  which  give  the  highest 
-■•r  -'rnricncy  when  primary  electrons  alone  are  con- 

‘•  'I  rD'  instance  the  parameters  are  so  adjusted  that 
‘  ii-iv  high  energy  particlea  hitting  the  central  or  Aat 

■  ■:  I 'lie  back  plate  which  is  at  the  highest  depressed 
•  ;.*■  Consequently  almost  all  the  back-scattered  elec- 

. ttiirr  reabsorbed  on  the  same  collector  that  they 

.  CO  to  another  collector  which  ia  one  stage  lower. 
'iMiiitl  losses  are  thua  minimized 
‘■K-  present  case,  the  trajectories  of  back-scattered 
••..•i-H  ;.re  traced  by  re-running  the  Herrmannsfeldt  code 
!  impacted  collectors  were  'cathodes’.  They  are 
‘“••I  ti:>t  outside  the  surface  of  the  electrode.  Figure 
.'.v«  r|i«.  trajectories  of  the  back-scattered  electrons  cor- 
■"l  iiiic  to  the  primaries  shown  in  Figure  5.  Certain 
assumptions  regarding  the  energy  and  direction 
ini'k-Mcattered  electrons  are  made,  in  order  to  obtain 

■  ’  r  .f  order  effects.  The  angle  of  'reffection'  is  made  equal 

■‘Hide  of  ’incidence’  with  respect  to  the  normal  to  the 
-•‘V’e.  {io|  The  energy  of  the  true  secondaries  is  assumed 

•  ••  -'ll)  fv.  The  energy  of  the  reAected  primaries  is  taken  to 
'  "  3  times  the  net  energy  of  impact  and  the  back-scatter 

‘^‘iit  is  taken  to  be  0.3  limes  the  primary  current,  these 
'  riHi  iM'jng  approximations  for  average  coefficients  in  the 
"  '•{  copper  as  the  material  for  the  electrodes.  [11]  In 
'"■•I  practice  the  angles  and  energies  would  have  a  natu- 

•  •;’r*‘.‘id.  The  back-scatter  coefficient  would  be  dependent 
*'  'pvin  the  surface  treatment  or  any  deposited  layer,  such 

-■riMin,  ns  well  as  the  angle  of  incidence. 

I  Ix’  trajectories  in  Figure  6  show  that  out  of  25  particles 
almost  all  the  se^daries  are  suppressed.  These 
tlx*  small  squiggles  at  the  points  of  impact  of  primaries, 
"v  are  returned  to  the  collector  with  the  next  highest 
"f*’*^d  potential;  and  one  returns  towards  the  entrance 
‘  Mie  collector  region.  Similarly  ali  the  reAected  primaries 
'Hlierted  on  the  same  electrode,  with  the  exception  of 
Mint  go  to  a  less  depressed  collector.  The  exceptional 


trajectories  happen  to  lie  at  tlie  edges  uf  ilie  energy  dis¬ 
tribution  curve,  and  have  less  wetghtage  liian  tlie  median 
particles.  Thus,  the  effect  of  these  bark-scattered  electrons 
in  reducing  the  estimated  collector  efficiency  is  minimized. 
With  the  present  level  of  approximations,  the  collector  ef- 
Aciency  is  estimated  to  lie  in  tlie  rage  of  70  to  80%. 

0.  CONCLUSIONS 

Results  of  a  proof-of-principle  experiment  have  been 
presented  for  a  double  cusp  system  for  energy  recovery  from 
spent  beam  in  a  large-orbit  gyrotron.  They  lend  support 
to  the  basic  approach  of  unwinding  the  beam  at  the  sec¬ 
ond  cusp.  A  depressed  collector  system  has  been  designed, 
including  trajectories  for  primary  as  well  as  back-scattered 
electrons.  Tins  system  has  demonstrated  energy  sorting, 
with  an  estimated  collector  efficiency  in  the  range  of  70  to 
80%. 
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Fig.  6.  Trajecloties  of  back-icattered  electron*  for 
parameters  and  points  of  impact  of  primaries  in  Figure  5. 
The  continuous  lines  are  reflected  primaries  and  the 
dotted  lines  are  true  secondaries. 
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The  propagation  of  high-power  ( 10-200  kW/cm^)  short-burst  (3-30  ns)  microwave  pulses  in 
the  atmosphere  has  been  studied  experimentally.  Microwave  power  from  a  large  orbit  gyrotron 
operating  at  9.6  GHz  is  focused  by  a  large-diameter  parabolic  reflector  into  a  test  cell.  The 
ambient  pressure  in  the  test  cell  was  varied  over  a  wide  range  and  the  microwave  power  density 
necessary  for  atmospheric  breakdown  has  been  determined  as  a  function  of  ambient  pressure 
and  pulse  duration.  Measurements  of  the  microwave  pulse  duration  before  and  after 
breakdown  have  been  obtained  to  determine  the  extent  to  which  microwave  energy  is  absorbed 
or  reflected  by  the  breakdown  plasma.  Results  are  compared  with  available  theory  and 
previously  reported  experiments. 


I.  INTRODUCTION 

Although  microwave  breakdown  of  the  atmosphere  has 
been  a  subject  of  study  for  many  years, considerable  new 
interest  has  recently  developed  in  the  field  arising  from  the 
availability  of  very  high-power  microwave  sources  devel¬ 
oped  for  directed-energy  systems  and  high-power  radar  ap¬ 
plications.  In  particular,  data  obtained  from  experiments  in¬ 
vestigating  microsecond-pulse  rf  breakdown  in  air  have 
supported  theoretical  predictions  that  very  high  microwave 
power  levels  could  be  propagated  in  the  atmosphere  if  the 
pulse  duration  could  be  made  sufficiently  short.  It  has  also 
been  proposed  that  at  low  ambient  pressures  a  very  high- 
power  microwave  pulse  might  accelerate  free  electrons  to 
energies  in  excess  of  that  associated  with  the  peak  in  the 
ionization  cross  section  for  air  in  a  distance  short  compared 
to  a  mean  free  path.  In  this  manner  even  higher  microwave 
powers  might  be  propagated  without  breakdown. 

Considerable  theoretical  and  experimental  work  has 
been  reported  in  this  area  in  recent  years. Theoretical 
results  reported  by  Woo  and  DeGroot, "  Yee  er  a/.,  and  Ali 
and  Coffey'^  have  indicated  that  when  breakdown  occurs, 
an  ionization  front  is  rapidly  formed  which  moves  toward 
the  microwave  source  and  consequently  decouples  the 
microwaves  from  the  initial  breakdown  region.  In  addition, 
the  amount  of  energy  that  can  be  transmitted  through  the 
medium  was  found  to  be  a  strong  function  of  the  initial  ener¬ 
gy  of  the  pulse  and  its  frequency,  pulse  shape,  and  duration. 

Several  experiments  have  been  reported  in  which  the 
propagation  of  short-burst  microwaves  through  air  has  been 
studied.  '*■”  Bollen  et  al.'*  have  reported  studies  of  microsec¬ 
ond-pulse  microwave  breakdown  in  nitrogen  by  focusing  the 
output  of  a  1 12-kW,  3S-GHz  gyrotron  into  a  test  cell.  Break¬ 
down  was  observed  at  microwave  power  densities  of  about 
30  kW/cm^  at  ambient  pressures  below  75  Torr.  Yee  el  al.  ’’ 
have  reported  studies  of  short-pulse  (10  ns)  microwave 
breakdown  of  air  in  a  rectangular  waveguide  powered  by  a 
1 5-MW  5-band  klystron  operating  at  2.856  GHz.  Results  of 
this  study  were  found  to  be  in  good  agreement  with  theory 
and  simulation.  Experimental  results  of  Gould  and  Ro¬ 
berts'^  and  Tetenbaum,  MacDonald,  and  Bandel'^  em¬ 
ployed  5-band  sources  with  repetition  rates  of  1200  and  20 


pulses  per  second,  respectively.  In  these  cases  the  test  cell 
was  a  resonant  cavity  fed  by  waveguide.  Finally,  recent  work 
by  Armstrong  et  al.'^  has  been  reported  in  which  high-power 
microwaves  from  an  5-band  relativistic  magnetron  have 
been  used  to  study  breakdown  thresholds.  In  these  experi¬ 
ments,  6-ns  microwave  pulses  were  focused  by  a  parabolic 
reflector  into  a  test  cell  and  the  resultant  breakdown  ob¬ 
served  using  optical  diagnostics. 

In  this  paper,  we  report  the  results  of  an  extensive  ex¬ 
perimental  study  of  short-burst  high-power  microwave 
propagation  through  air.  Data  have  been  obtained  for  trans¬ 
verse  electromagnetic  (TEM)  microwave  pulse  propagation 
in  a  test  cell  large  enough  to  ensure  that  surface  breakdown 
does  not  occur  to  any  significant  extent.  Furthermore,  data 
have  been  obtained  over  a  wide  range  of  microwave  pulse 
power  densities  and  pulse  durations  and  for  an  equally  wide 
range  of  ambient  test-cell  pressures.  Section  II  contains  a 
review  of  available  theory  relating  to  these  experiments.  Sec¬ 
tion  III  discusses  the  experimental  apparatus  used  for  the 
study,  and  results  are  presented  in  Sec.  IV.  Conclusions  are 
drawn  in  Sec.  V. 

II.  THEORETICAL  DISCUSSION 

The  theory  of  microwave  breakdown  of  a  neutral  media 
has  been  extensively  treated  in  early  work  by  MacDonald.^ 
It  is  assumed  that  breakdown  occurs  when  the  electric  field 
of  the  wave  accelerates  free  electrons  in  the  media  to  energies 
sufficient  to  ionize  neutral  species  and  initiate  an  avalanche 
process  that  produces  a  rapid  increase  in  the  plasma  density. 
The  continuity  equation  for  electrons  in  the  media  can  be 
written 

~  =  Vin —  v,n +  7^{Dn),  (1) 

dt 

where  n  is  the  electron  density,  v,  is  the  ionization  rate,  v.  is 
the  attachment  rate,  and  D  is  the  electron  diffusion  coeffi¬ 
cient.  The  diffusion  term  can  be  approximated  as 
V'(Z7n)  afl(  —  D/S}),  where  A  is  the  characteristic  diffu¬ 
sion  length  and  depends  upon  the  system  geometry.  Using 
this  approximation,  the  continuity  equation  has  a  solution 
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where  rtg  is  the  initial  free-electron  density.  The  diffusion 
coefficient  D  was  found^  to  follow  the  empirical  relation 

Dp  =  [29 +  0.9(E/p)]X\0*,  (3> 

where p  is  the  ambient  media  pressure  and  E  is  the  amplitude 
of  the  wave  electric  field.  MacDonald  noted  that  in  cw  mi¬ 
crowave  breakdown  phenomena  steady  state  is  reached  and 

v,-v,  =Z)/A^  (4) 

He  then  matched  these  results  to  data  obtained  from  cw  mi¬ 
crowave  breakdown  experiments  conducted  in  a  cavity.  In 
this  manner  he  obtained  estimates  of  the  net  ionization  rate 
V,  —  as  a  function  of  ambient  media  pressure  for  the  given 
wavelength  and  electric  field  amplitude. 

For  pulsed  microwave  breakdown,  Eq.  (4)  does  not  ap¬ 
ply  and  the  electron  density  grows  at  a  rate  given  by 

v  =  v,-v,-D/A\  (5) 

A  reasonable  estimate  of  the  time  taken  to  achieve  break¬ 
down  will  be  obtained  if  it  is  assumed  that  breakdown  occurs 
as  the  plasma  density  approaches  the  critical  plasma  density, 

/Ip  =  10'V/i*cm~^  (6) 

The  pulse  duration  r  for  breakdown  to  occur  can  then  be 
determined  approximately  from 

zip  =  /loe"  (7) 

Experimentally,  the  initial  free-electron  density  is  not  usual¬ 
ly  known  to  any  degree  of  accuracy,  but  the  exponential  na¬ 
ture  of  the  rise  in  the  electron  density  with  time  tends  to 
minimize  the  effect  of  errors  in  the  number  assumed  for  /Iq  . 
In  addition,  studies  by  Rappaport,  Latham,  and  Striffler'* 
have  indicated  that  wave  propagation  can  be  significantly 
affected  by  plasma  densities  considerably  below  the  critical 
plasma  density,  although  the  exponential  nature  of  the  den¬ 
sity  increase  with  time  minimizes  errors  resulting  from  this 
assumption  as  well. 

In  a  similar  approach,  Ali  and  Coffey'^  theoretically 
studied  breakdown  by  treating  electrons  with  a  Maxwellian 
velocity  distribution.  In  this  work  the  average  electron  ener¬ 
gy  was  assumed  to  be  linearly  dependent  on  E  /p.  They  as¬ 
sumed  that  the  ionization  in  air  is  primarily  from  the  oxygen 


molecule  and  used  known  values  for  the  attachment  fre¬ 
quency  and  estimates  for  the  diffusion  scale  length  A  to  de¬ 
termine  the  net  ionization  rate  as  a  function  of  ambient  pres¬ 
sure  and  pulse  duration. 

In  addition  to  this  analytical  work,  significant  numeri¬ 
cal  simulation  has  been  conducted  of  short-pulse  microwave 
propagation  in  the  atmosphere, results  in  reasonable 
agreement  with  experiment  where  data  are  available. 

ill.  EXPERIMENTAL  APPARATUS 

High-power  microwave  pulses  used  for  the  propagation 
studies  were  generated  by  a  large-orbit  gyrotron’’  "  operat¬ 
ing  at  10  GHz.  This  device,  powered  by  a  pulse  line  accelera¬ 
tor  [2  MeV,  10-20  kA,  30  ns  full  width  at  half-maximum 
(FWHM)],  produces  high-power  microwave  radiation 
( 100-500  MW)  by  the  resonant  interaction  of  the  modes  of 
an  axis  encircling  rotating  electron  beam  with  the  modes  of  a 
magnetron-type  multiresonator  circuit.  The  rotating  elec¬ 
tron  beam  is  produced  by  passing  a  cylindrical,  nonrotating 
beam  from  a  field  emission  diode  through  a  narrow  magnetic 
cusp,  as  shown  in  Fig.  1.  The  cusp  acts  to  convert  the  up¬ 
stream  electron  axial  velocity  (i/^, )  into  a  downstream  azi¬ 
muthal  (1/02 )  and  axial  velocity  ^  manner  given 

approximately  by 

+  >^82  =  vh  +  (8) 

Here  =  eB  /my  is  the  relativistic  electron  cyclotron  fre¬ 
quency,  where  e  is  the  electronic  charge,  m  is  the  electron 
rest  mass,  B  is  the  axial  magnetic  field,  and 

y=  (1  _y2/c2)-''2  (9) 

is  the  relativistic  mass  factor  for  the  beam  electrons  (c  is  the 
speed  of  light  in  vacuum ) .  The  radius  is  the  mean  radius  of 
the  rotating  electron  beam,  which  for  a  symmetric  cusp  is 
approximately  equal  to  the  cathode  radius.  Thus  for  a  given 
applied  axial  magnetic  field,  only  electrons  with  energies 
greater  than  a  value  given  by 

=  +  (10) 

can  pass  through  the  cusp  into  the  downstream  drift  region. 

The  above  result  allows  for  the  use  of  this  facility  for  the 


"Magnetron"  Outer  Conducting 


FIG.  1 .  Large  orbit  gyrotron  experimental  coo' 
figuration. 
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production  of  very  short  microwave  pulses.  Since  the  diode 
voltage  waveform  slowly  rises  to  a  maximum  value  and  then 
slowly  decays,  the  duration  of  the  downstream  electron 
beam  pulse  can  be  varied  within  certain  limits  by  varying  the 
applied  axial  magnetic  field  strength.  In  this  manner  ail  elec¬ 
trons  with  energies  below  that  value  required  for  transmis¬ 
sion  through  the  cusp  are  reflected  and  the  downstream 
beam  pulse  can  be  considerably  shorter  in  duration  than  the 
diode  pulse,  as  shown  in  Fig.  2.  Although  optimum  micro- 
wave  production  cannot  be  achieved  over  as  wide  a  variation 
of  magnetic  field  as  might  be  desired,  data  have  been  ob¬ 
tained  for  microwave  pulses  of  varying  power  in  the  range  3- 
30  ns  (FWHM). 

The  experimental  configuration  used  for  the  propaga¬ 
tion  experiments  is  shown  in  Fig.  3.  Microwaves  from  the 
large-orbit  gyrotron  are  focused  by  a  large  diameter  (150 
cm)  parabolic  reflector  with  a  75-cm  focal  length  into  a  30- 
cm-diam  cylindrical  acrylic  vacuum  chamber  which  served 
as  the  test  cell.  The  test  cell  was  fitted  with  a  high-accuracy 
vacuum  gauge  and  a  calibrated  leak  valve  that  allowed  main¬ 
tenance  of  ambient  pressures  from  atmospheric  down  to  0. 1 
Torr.  Measurements  of  the  microwave  power  density  along 
the  dish  axis  within  the  test  cell  are  shown  in  Fig.  4.  The 
points  plotted  are  the  averages  of  four  shots  taken  at  each 


FIG.  3.  Configuration  used  for  microwave  propagation  studies. 


FIG.  4.  Microwave  power  density  vs  axial  distance  from  the  crater  of  the 
test  cell.  The  test<ell  diameter  is  30  cm. 


position.  In  this  way  the  effects  of  shot-to-shot  variation  in 
the  microwave  output  power  of  the  device  can  be  minimiypH 
Although  the  focusing  of  the  microwaves  into  the  test  cell 
was  not  as  effective  as  had  been  expected  from  low-power 
measurements,  it  was  sufficient  to  ensure  that  breakdown, 
when  observed,  did  not  initiate  on  the  test-cell  walls.  Surface 
breakdown,  which  is  generally  easier  to  initiate  than  bulk 
media  breakdown,  must  be  avoided  if  true  media  breakdown 
limits  are  to  be  established.  The  superposition  of  the  incident 
wave  from  the  source  with  the  reflected  wave  from  the  para¬ 
bolic  reflector  results  in  a  small-amplitude  standing  wave  in 
the  cell,  although  measurements  indicate  that  the  incident 
microwave  power  is  more  than  a  factor  of  10  lower  than  the 
reflected  microwave  power  at  the  focus.  No  seed  ionkation 
source  was  used  for  these  studies. 


FIG.  3.  Schematic  of  microwave  diagnostic  setup  used  in  the  propagation 
experiments. 
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FIG.  6.  Time-inlegrated  photographs  of  microwave-induced  breakdown  in 
the  lest  cell,  (a)  4.8,  (b)  9.S,  (c)  20.  and  (d)  28  Torr. 

Open-ended  .f-band  (8-12  GHz)  waveguide  sections 
were  used  to  detect  the  microwave  signals  before  and  after 
the  focal  point  in  the  test  cell,  as  shown  in  Fig.  S.  The  micro- 
wave  signals  were  passed  through  calibrated  fixed  attenua¬ 
tors  to  calibrated  detectors.  All  components  were  calibrated 
inhouse  using  standard  techniques.  The  detectors  chosen 
were  selected  on  the  basis  of  their  measured  transient  pulse 
response  time  (rise  time  ~  1  ns).  Detector  signals  were  fed 
directly  into  Tektronix  7104  fast  oscilloscopes. 

IV.  EXPERIMENTAL  RESULTS 

In  previous  work,'^  the  microwave  pulse  duration  that 
could  be  propagated  without  breakdown  was  determined  by 
simply  me;;;iuring  the  duration  of  the  post-breakdown  mi¬ 
crowave  pulse.  The  experiments  reported  here  were  de¬ 
signed  so  that  breakdown  could  be  detected  by  two  indepen¬ 
dent  criteria.  The  pulse  durations  before  and  after  the  focal 
point  in  the  test  cell  were  compared  to  determine  if  micro- 
wave  absorption  and  reflection  by  breakdown  plasma  had 


(b) 


FIG.  7.  Typical  microwave  detector  wavefoms,  (a)  incident  detector,  and 
(b)  transmitted  detector.  Waveforms  correspond  to  case  (a)  in  Fig.  6. 

shortened  the  duration  of  the  post-focus  microwave  pulse 
significantly.  In  addition,  time-integrated  photographs  of 
the  visible  light  produced  by  ionization  processes  in  the  test 
cell  were  obtained  using  a  still  camera.  A  typical  set  of  pho¬ 
tographs  of  breakdown  in  the  test  cell  is  shown  in  Fig.  6,  and 
they  clearly  indicate  that  surface  breakdown  at  the  test-cell 
wall  is  not  a  significant  factor  in  the  measurements.  For  a 
fixed  microwave  power  density  and  pulse  duration,  the 
breakdown  region  becomes  progressively  smaller  as  the  am¬ 
bient  test-cell  pressure  is  increased.  At  the  higher  pressures 
evidence  of  the  small-amplitude  standing  wave  in  the  test 
cell  can  be  clearly  seen  in  the  photographs.  The  observed 
wavelength  of  the  breakdown  pattern  is  fully  consistent  with 


FIG.  8  Microwave  propagation  data  plotted  using  the  visible  light  criterion 
for  breakdown. 
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FIG.  9.  Microwave  propagation  data  plotted  using  the  pulse  width  criterion  FIG.  10.  Comparison  of  results  of  present  study  with  available  theory  and 

for  breakdown.  experiment. 


the  operating  frequency  (9.6  GHz)  of  the  large-orbit  gyro- 
tron. 

A  typical  set  of  oscilloscope  waveforms  for  the  prefocus 
and  post-focus  detector  systems  are  shown  in  Fig.  7.  The 
examples  shown  are  for  a  case  where  breakdown  resulted  in 
considerable  shortening  of  the  post-focus  microwave  pulse 
duration.  Data  from  all  experiments  are  plotted  in  the  famil¬ 
iar  E  /p  vs  pr  format  in  Figs.  8  and  9.  The  data  points  in  the 
two  plots  are  the  same,  but  in  Fig.  8  the  visible  light  criteria 
for  breakdown  was  used  to  indicate  on  which  shots  break¬ 
down  occurred,  and  in  Fig.  9  breakdown  was  defined  by  a 
measured  post-focus  pulse  duration  of  80%  or  less  of  the  pre¬ 
focus  pulse  duration. 

V.  CONCLUSIONS 

From  the  data  it  is  possible  to  draw  a  line  on  the  £  /p  vs 
pr  plot  above  which  breakdown  always  occurs  and  a  second 
line  below  which  breakdown  never  occurs.  These  results 
have  been  compared  with  data  from  experiments  by  Yee  et 
Gould  and  Roberts,'^  Tetenbaum  et  al.,'*  and  Fel- 
senthal  and  Proud'”  in  Fig.  10.  Also  shown  are  the  predic¬ 
tions  of  theory  by  MacDonald’  and  Ali  and  Coffey”  for 
wavelengths  of  3  and  10  cm. 

The  data  obtained  in  these  experiments  are  in  reasona¬ 
ble  agreement  with  data  obtained  by  Yee  et  for  break¬ 
down  inside  a  waveguide  at  2.85  GHz.  The  results  are  also  in 
reasonable  agreement  with  the  theory  of  MacDonald.’  Both 
the  present  data  and  the  data  of  Yee  er  al.  fall  somewhat 
below  the  theoretical  predictions  of  Ali  and  Coffey,  although 
the  qualitative  agreement  between  theory  and  experiment  is 
very  good.  It  is  interesting  to  note  that  no  seed  ionization 
source  was  employed  in  these  experiments,  in  contrast  to  the 
experiments  reported  by  Yee  et  a!.''  and  Gould  and  Ro¬ 
berts.”  Usually,  seed  ionization  sources  are  used  to  provide  a 
reliable  source  of  small  numbers  of  free  electrons  to  impove 
reproducibility.  In  the  experiments  reported  here,  reproduc¬ 
ibility  was  not  a  major  problem,  although  it  is  likely  that 


X  rays  produced  by  energetic  electrons  striking  the  drift  tube 
walls  produced  sufficient  free  electrons  to  initiate  the  break¬ 
down  process. 

Future  experiments  are  planned  to  extend  the  range  of 
data  obtained  to  even  higher  microwave  power  densities  and 
shorter  pulse  durations.  In  addition,  the  propagation  of 
high-power  microwave  pulses  through  preformed  plasmas  is 
also  under  study. 
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Propagation  and  Breakdown  in  Gases  by  Short  Burst,  High-Power  Microwaves* 
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20742 


We  are  studying  propagation  of  pulsed  high-power  microwaves  in  the  atmosphere.  Indi¬ 
vidual  processes  in  our  problem  include  reflection,  absorption,  and  transmission  of  normadly 
incident  electromagnetic  radiation  onto  a  finite  width  resistive  plasma  and  ionization  of  vari¬ 
ous  gases  by  high-power  RF.  Wc  examine  methods  to  solve  for  propagation  characteristics  of 
a  time- varying  inhomogeneous  plasma  where  the  plasma  density  profile  evolution  depends  on 
the  amplitude  of  the  transmitted  electric  field.  The  growth  of  a  plasma  in  a  one-dimensional 
gas  cell  using  a  low  ionization  rate,  long-pulse  model  is  found.  The  validity  of  this  model 
is  verified  by  fast  Fourier  transform  methods.  Also  discussed  is  a  method  for  solving  for 
propagation  characteristics  of  a  time-varying  plzisma  half-space  in  the  critical  region  where 
geometrical  optics  techniques  break-down. 

*This  work  is  supported  by  AFOSR. 
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STATEMENT  OF  PROBLEM 


PURPOSE  -  TO  UNDERSTAND  THE  PROPAGATION  PROPERTIES  OF  HIGH  POWERy 
SHORT  PULSE  MICROWAVES  THROUGH  GASES. 


CONFIGURATION  AND  PARAMETERS  OF  EXPERIMENT  AT  MARYLAND 
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PROPAGATION  OF  RF  IN  GASES 
AND  PLASMAS 


Problem  Overview 


The  total  problem  of  propagation  and  breakdown  in  gases  may  be  di¬ 
vided  up  into  the  individual  processes  of  ionization  of  gases  by  RF  and 
propagation  of  RF  through  a  time-varying  inhomogeneous  plasma.  We 
take  the  view  here  that  if  we  can  solve  the  propagation  problem  we  Ccin 
solve  the  total  problem  by  demanding  self-consistency  of  the  time-evolving 
inhomogeneous  plasma  density  profile  with  the  calculated  ionizing  electric 
field  propagating  through  the  plasma.  The  connection  between  the  electric 
field  and  the  plasma  density  is  made  through  ionization  and  attachment 
rate  equations  for  air.^  The  motivation  for  the  overall  approach  lies  in  the 
fact  that  the  total  problem  of  ionization  and  propagation  combined  is  a 
non-linear  problem.  By  considering  the  propagation  problem  and  ioniza¬ 
tion  problems  separately  we  first  solve  a  linear  wave  equation  with  varying 
coefficients  and  then  afterwards  impose  self-consistency  on  the  density  evo¬ 
lution.  This  might  be  done  by  trial  and  error,  a  relaxation  method  or  an 
iterative  technique. 

In  the  first  portion  of  this  poster  we  consider  a  parameter  regime  where 
the  ionization  rate  of  the  plasma  is  so  low  that  plasma  ionization  and  RF 
propagation  are  on  different  time  scales.  In  this  case  methods  for  solving 
for  the  electric  field  in  an  inhomogeneous  time-independent  plasma  density 
are  used  and  self-consistcncy  with  the  ionizing  electric  field  is  imposed. 


In  the  second  portion  of  this  poster  the  propagation  problem  for  a  pa¬ 
rameter  regime  where  the  density  may  not  be  considered  static  is  discussed. 
Here  the  ionization  rate  is  still  held  very  low  compared  with  the  RF  fre¬ 
quency  and  so  the  term  medium  ionization  rate  model  seems  appropriate. 

In  future  work  we  may  study  a  parameter  regime  where  the  ionization 
rate  is  of  the  order  of  the  RF  frequency.  This  would  be  called  a  high 
ionization  rate. 

^W.  Woo  and  J.S.  Degroot,  Phys.  Fluids  27,  475  (19S4). 


Low  Ionization  Rate  Model  for  Microwave  Breakdown 

of  Air 


We  consider  the  problem  of  propagation  of  a  normally  incident  plane 
wave  through  a  1-dimensional  gas  cell.  To  be  found  are  transmitted  and 
reflected  field  amplitudes  as  functions  of  time.  Assuming  that  the  ioniza¬ 
tion  rate  is  much  smaller  than  the  carrier  frequency  Uc  and  that  electron 
momentum  change  due  to  ionization  is  negligible  compared  to  momentum 
transfer  due  to  elastic  scattering,  we  model  our  plasma  as  follows: 
Maxwell: 


V^E 


1  d~E  dJ 

c2  dt 


where  J  =  qnv 
Momentum  balance: 


Electron  continuity: 


Charge  neutrality: 


dv 

m—  =  gE  -  mvcv 


dn 

dt 


=  -  i'a)n 


^  —  ^e{i)  —  ^  ^  ^a^a(ion8)(^) 


where  =  1  for  singly  ionized  molecules,  2  for  doubly  ionized  molecules 


etc. 


Ionization,  attachment,  and  elastic  scattering  rates  for  air  from  Woo  and 
Degroot; 


Vi  =  Ptoff  5.14  X  10^^  exp(— TSof 


-.44' 


32  <  q;  <  100 


=  Ptoff  7.6  X  10~^a2(c^  -f  218) 
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PToFFv/l+a;2/t/2 

In  this  poster  Va  is  taken  to  be  zero  when  60  <  a  <  100.  Solution  of  the 
electron  continuity  equation  yields: 

n(t)  =  n(0)e<‘'‘-‘'->' 

The  momentum  balance  equation  is  solved  with  an  integrating  factor  to 
yield: 

u  =  ~  /  dr E(t) 

m  J 

So  the  1-dimensional  wave  equation  becomes 
d^E 


dz^ 


where  ujp{t)  =  n(t)q^/m€Q  and  again  the  ionization  rate  is  tsiken  to  be 
much  smaller  than  the  carrier  frequency.  If  we  now  further  assume  that  the 
electromagnetic  response  time  =  Tem  of  the  gas  cell  is  much  shorter  than 
1/vi  we  may  take  solutions  to  the  wave  equation  of  the  form 


for  a  spatially  uniform  plasma.  The  appropriate  dispersion  relation  is 

1  -  tUc/ujc 
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PROPAGATION  OF  NORMAL  INCIDENT  PLANE  WAVE 
THROUGH  A  FINITE  THICKNESS  REGION  OF  IONIZED  AIR 
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By  matching  boundary  conditions  we  find 

Et  ^  -4npe>*^»^ 

(1  —  rip)^  ~  (1  + 

~  (1  -  np)2  -  (1  +  np)2e2>fc.«i 


where 


Discussion  and  Results 


The  parameters  chosen  for  numericail  computation  of  results  in  this 
problem  axe  incident  wave  amplitude  =  1100  V/cm,  carrier  frequency  fc 
=  10  Ghz,  gas  cell  length  =  4  cm,  and  gas  cell  pressure  =  5  Torr.  In  the 
first  set  of  plots  a  uniform  density  model  is  used,  the  electron  density  in 
the  gas  cell  is  assumed  to  be  spatially  uniform  but  increasing  with  time. 
At  each  time  step  of  the  calculation  the  transmitted  field,  reflected  field 
and  the  field  at  the  center  of  the  gas  cell  are  found.  The  field  at  the  center 
of  the  gas  cell  is  used  to  calculate  the  ionization  and  attachment  rates  for 
the  gas  which  in  turn  are  used  to  find  the  electron  density  for  the  next 
time  step.  We  see  that  the  transmitted  wave  falls  sheirply  in  amplitude  as 
the  density  approaches  the  critical  density  ric  =  uj\m€Q/q^  at  t  «  250  ns. 
The  initial  value  of  n/ric  —  10~®  corresponds  to  an  initial  electron  density 
of  10^  electrons/cm^.  Note  that  Vi  10®  zind  Uc  10^°  so  i/,-  <C  as 
required.  Also  Tem  ~  .3  ns  as  will  be  discussed  momentarily  so  l/i/j  ^  Tem 
as  required.  Duration  of  each  time  step  is  1  ns. 

In  the  next  two  plots  a  twenty  step  model  is  used  and  the  gas  cell  is 
broken  up  into  twenty  4  cm/20  =  .2  cm  regions.  The  density  in  each  region 
is  calculated  separately  at  the  beginnii.g  of  each  time  step.  A  set  of  42 
equations  in  42  unknown  variables  is  solved  numerically  to  find  forward  and 
reverse  field  amplitudes  in  each  region  as  well  as  reflected  and  transmitted 
wave  amplitudes  outside  the  gas  cell.  The  electric  field  amplitude  at  the 


center  of  each  region  is  used  to  calculate  the  ionization  rate  for  that  region 
at  each  time  step.  Shown  is  the  evolution  of  the  density  profile  and  the 
evolution  of  the  electric  field  amplitude  across  the  gas  cell  with  time.  Each 
time  step  lasts  1  ns.  Note  that  in  these  plots  we  may  think  of  the  plane 
wave  driving  the  gas  cell  as  being  incident  from  the  left.  The  electric  field 
amplitude  decreases  across  the  gas  cell  because  of  elastic  collisions.  This 
in  turn  causes  a  lower  ionization  rate  on  the  right  side  of  the  gas  cell  and 
plasma  buildup  is  greatest  on  the  left  side. 

The  final  set  of  plots  allows  us  to  evaluate  •  Here  the  Fourier  spec¬ 
trum  of  an  input  pulse  is  multiplied  by  the  frequency  response  of  a  uniform 
in  space,  constant  in  time,  density  plasma  of  length  4  cm,  P  =  5  Torr  and 
»^/”carrier  =  -5  .  ncarrier  is  the  critical  plasma  density  at  /  =  /carrier  =  10 
Ghz.  The  output  pulse  in  frequency  and  time  domains  is  shown.  Since  the 
input  pulse  starts  at  t  =  3  ns  and  the  output  pulse  reaches  its  steady  state 
value  at  i  «  3.3  ns,  ^  .3  ns.  A  similiar  value  of  is  seen  for  the  case 
of  a  4  cm,  5  Torr,  n/ncarrier  =  1.0  plasma. 


Problem  Parameters 


Carrier  frequency  =  10  Ghz 

Gas  cell  length  =  4  cm 

Gas  cell  pressure  =  5  Torr 

Initial  electron  density  =  1000  electrons/cm^ 

For  uniform  density  model  and  twenty  step  model  incident 
plane  wave  amplitude  =  1100  V/cm 


Relevant  Time  Scales 

l/i'i  =  1.5  X  10"®  seconds  for  Eras  =  1100  V/cm 
'^em  3  X  10~^®  seconds 
1/^c  =  1-6  X  10~^^  seconds 
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Medium  Ionization  Rate  Model  for  Microwave  Propagation 

through  Plasma 


In  the  low  ionization  rate  model  we  assumed  that  the  ionization  rate 
Ui  -C  l/'Tem  where  Tem  is  the  electromagnetic  response  time  of  the  system 
under  consideration.  In  the  medium  ionization  rate  model  we  remove  this 
restriction  and  require  only  that  i/,-  -C  uJc-  The  geometrical  optics  method^ 
permits  solution  of  the  wave  amplitude  in  a  time-varying  inhomogeneous 
plasma.  This  method  converges  fastest  as  A  0.  In  a  static  homoge¬ 
neous  collisionless  plasma  where  Up  «  Uc  (at  critical  density)  A  — >•  oo  and 
geometrical  optics  breaks  down.  Thus  another  method  is  needed. 


Motivation  for  Power  Series  Method 

In  a  static  collisionless  1-dimensional  plasma  the  wave  equation 

d^E  1  d^E  ujI  „ 

dz^  C2  dt^  C2 

becomes  in  a  critical  density  region 

d^E 


az2 


«  0 


This  suggests  a  perturbation  expansion  about  the  critical  density 

E  —  Eq  -j-  El  -j-  E2  *4"  •  •  • 


d^Eo 


where 


and 


Eo{z,t)  ^  Ait)z  +  B{t) 

Consideration  of  the  first  order  equation 

1  d^Eo 

dz^  <?  dfi  ^  c2  “ 

shows  US  that  to  first  order  the  solution  for  £"  is  a  cubic  in  z  with  coefficients 
that  are  functions  of  t.  To  second  order  the  solution  for  £*  is  a  fifth  degree 
polynomial  in  z  with  coefficients  that  are  functions  of  t.  Thus  expanding 
to  all  orders  we  have  a  solution  for  E  which  is  a  power  series  in  z  with 
coefficients  that  are  functions  of  t.  Convergence  is  expected  to  be  best  in 
regions  where  Wp  «  cuc,  i-e.  the  region  where  geometrical  optics  fails. 

Bernstein,  Phys.  Fluids  IS,  320  (1975). 


Propagation  of  Normal  Incident  Plane  Wave 
Onto  Time  Varying  Plasma  Half-Space 


free  space 

time-varying  density  homogeneous  plasma 

incident  plane  wave 

transmitted  plane  wave 

nnrryyrrrrn^ 

tn'Trrrrrrx^ 

Ei 

Et 

reflected  plane  wave 

^>TrrvYVTVVT> 

Er 

plugging  into  wave  and  momentum  balance  equations  and  equating  like 
powers  of  z  gives  us  a  pair  of  recursion  relations  for  the  an{t)  and  bn{t) 

an+2(0(^  +  l)(n  +  2)  =  ~ 

dhn  , 

— —  — •  dn  ^cb-n 
at 

By  truncating  the  series  after  iV  -f- 1  terms  where  N  is  the  degree  of  the 
resulting  polynomial  a  set  of  2(Ar  +  l)  coupled  second  order  O.D.E.’s  result. 
These  are  easily  solved  numerically  by  a  Runge-Kutta  code  provided  suit¬ 
able  initial  conditions  Un(0),<2|i(0),6n(0),6Jj(0)  are  found.  An  e.xponentially 
increasing  density  is  of  interest: 

ne(z,t)  =  ne(O)0(z)  exp(i/ii  0(i)) 
where  the  0  function  is  defined 


^(77)  =  1  77  >  0 

0(tj)  =  0  77  <  0 

For  this  case  the  a„(0)  are  easily  found  since  the  plasma  is  in  steady  state 
condition  for  ^  <  0. 


£{2,  (  =  0)  =  ^  a„(0)»"  = 

n=0 


{-ikz)" 
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a„(0)  =  Et 
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where 


Et 


=  Ei 


2kQ 
k  kQ 


and 


Find  6„(0) 


^v{2,t  =  0)  = 

^  n=0 


In  steady  state  (t  <  0) 


iUrTTl  „ 

—^v  -  E  - 

q 


Uctn 
- V 
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(iojc  +  i^c)  — V  =  E 
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Therefore 


oo 


oo 


E'’->(o>'’  =  E 


Qn(0) 
iUc  +  Uc 


The  6^,(0)  are  given  by 


^-nCO) 


Qn(0) 
iUc  +  Uc 


KW  =  an(0)  -  l/e6n(0) 

Find  the  a'„(0)  from  the  following  boundary  conditions  for  the  electric  field 


dE\ 
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Write  the  electric  field  for  2  =  O'  as  the  sum  of  an  incident  term  and  a 


reflected  term  Fr 


E{z,t)  =  EiC  +  iv(^ -f  ct) 


2  <  0 


=  -i  — +  F^ct) 

9^  J  ,=o- 

+  cF^ct) 

9t  J  „o- 

combining  and  using  the  boundary  condition  produces 


=  2iWe£ie’'“'‘ 


/  zs=0+  /  z=0+ 


Introduce  the  power  series 


£'=^an(0"”  z>0 


and  we  have 


This  relationship  produces  ao(0)  immediately  from  ai(0) 


ao(0)  =  2iuJcEi  +  cax(O) 


repeated  differentiation  of  equation  *  and  liberal  use  of  the  recursion  re¬ 
lations  for  the  Gn  and  bn  permit  finding  all  of  the  n'i(O).  This  repeated 
differentiation  process  also  allows  us  to  find  the  a'/i(0),  a'//(0)  etc. 


The  above  method  finds  the  an(0),6n(0),a^(0),6'„(0)  which  permits  a 
Runge-Kutta  code  to  find  the  electric  field  as  a  function  of  2  for  z  >  0  at 
any  time  step.  Plots  shown  of  E  versus  2  are  made  from  this  method. 

Alternatively  if  we  wish  only  to  find  the  electric  field  for  2  =  O''’  (and 
from  the  boundary  conditions  the  electric  field  for  2  =  0~)  as  a  function  of 
time  we  may  use  the  derivatives  of  ao(t) 


“  a«t« 


£(z  =  0^^)  =  E^ 


t=0 


l\ 


Plots  shown  of  E  versus  t  are  made  from  this  method. 


Parameters  Used  in  Plotting 


E  versus  z  plot 

1.  n/nc)t=o  =  .5 

2.  Ui  —  .IcJc 

3.  Uc  =  0 

4.  Runge-Kutta  time  step  =  .S/c^c 

5.  Output  is  after  1  time  step  — *•  t  =  .b/uc 

6.  n/ric  =  n/nc)t=o  exp(iyit)  =  .5exp(.05)  =  .526 

7.  Degree  of  approximating  polynomial  =  20. 

E  versus  t  plot 

1.  n/nc)t=o  —  .5 

2.  Vi  —  Auc 

3.  i/c  =  0 

4.  Curve  labeled  actual  is  a  plot  of  ao(f)  with  32  terms  kept  in  the  taylor 
series. 

5.  Curve  labeled  static  is  a  plot  of 


TYPICAL  E  VERSUS  Z  PLOT 


Future  Work 


1.  Study  convergence  properties  of  power  series  solution  for  the  electric 
field  in  a  time- varying  plasma. 

2.  Investigate  the  physics  of  propagation  of  RF  through  a  time- varying 
plasma  which  has  density  increasing  beyond  critial  density. 

3.  Consider  alternative  expansions  for  the  electric  field  in  a  time- varying 
plasma  and  their  convergence  properties. 

4.  Impose  self-consistency  on  a  plasma  density  profile  with  an  ionizing 
propagating  electric  field  for  a  plasma  half-space  and  other  geometries. 

5.  Study  high  ionization  rate  breakdown  where  the  ionization  rate  is  of 
the  order  of  the  RF  carrier  frequency. 


Shielded-Source,  Short-Pulse  Microwave  Propagation  Experiments 


W.  W.  Destler,  Z.  Segalov,  and  J.  Rodgers 
Electrical  Engineering  Department 
University  of  Maryland,  College  Park  MD  20742 

Abstract 

The  propagation  of  short-pulse,  high  power  microwaves  through  a  neutral  gas 
has  been  studied  and  breakdown  limits  determined  in  an  experimental  configuration 
in  which  the  neutral  gas  test  cell  was  carefully  shielded  from  x-rays  emanating  from 
the  high  power  microwave  tube.  In  this  manner  the  possible  x-ray  generation  of 
free  electrons  in  the  test  cell  was  minimized.  Residts  are  compau-ed  with  previous 
experiments  in  which  the  test  cell  was  not  shielded  from  the  source. 


1 


I.  Introduction 

Renewed  interest  in  the  propagation  of  high  power  microwave  pulses  in  the 
atmosphere^”®  has  resulted  from  the  availability  of  very  high  power  microwave 
sources  with  pulse  durations  in  the  range  2-1000  ns/  It  has  been  suggested  that 
extremely  high  levels  of  microwave  power  can  be  propagated  in  the  atmosphere  if 
the  microwave  pulse  duration  is  significantly  shorter  thzm  the  time  required  to  build 
up  electron  densities  sufficient  to  absorb  or  reflect  significant  microwave  energy. 

In  previous  work®,  our  group  has  reported  studies  of  the  propagation  of  high 
power  (o0-300  kW/cm^),  short-pulse  (2-35  ns)  microwaves  through  an  air  filled  test 
cell.  In  these  experiments,  microwaves  from  a  high  power  Large  Orbit  Gyrotron  were 
focused  into  a  test  cell  using  a  large  diameter  parabolic  reflector,  and  microwave 
breakdown  thresholds  were  determined  over  a  wide  range  of  test  cell  pressure  and  mi¬ 
crowave  power  density  and  pulse  duration.  These  free  space  microwave  breakdown 
thresholds  were  found  to  be  comparable  to  previous  data  obtained  on  breakdown 
limits  in  waveguides,  and  shorter  pulse  durations  did  allow  higher  microwave  power 
densities  to  propagate  without  breakdown,  as  expected. 

It  has  been  suggested®,  however,  that  x-radiation  produced  by  the  energetic 
electrons  in  the  Large  Orbit  Gyrotron  could  be  responsible  for  generating  a  num¬ 
ber  of  free  electrons  in  the  test  cell,  zmd  that  breakdown  might  therefore  have 
been  observed  at  artificially  low  values  of  microwave  power  density.  In  the  exper¬ 
iments  reported  herein,  the  test  ceil  was  carefully  screened  from  the  Large  Orbit 
Gyrotron  x-radiation  by  a  concrete  and  lead  shield  wall,  and  breakdown  limits  were 
determined  in  the  manner  reported  previovisly.  Thus  a  reawly  comparison  between 
breakdown  limits  with  and  without  source  x-ray  shielding  is  now  possible. 

In  this  paper,  related  theoretical  considerations  are  discussed  in  section  11.  The 
experiments  conducted  are  detailed  in  section  III,  and  results  are  compared  with 
previously  obtained  unshielded  source  breakdown  limits  in  section  IV. 

II.  Theoretical  Discussion 

The  theory  of  microwave  breakdown  of  a  neutral  media  is  treated  extensively  in 
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early  work  by  MacDonald’®.  T’.iis  treatment  assumes  that  breedcdown  occurs  when 
the  electric  field  of  the  microwaves  accelerates  any  free  electrons  in  the  media  to  en¬ 
ergies  sufficient  to  cause  ionizir.g  collisions  with  neutrals.  The  resulting  production 
of  additional  electrons  drives  an  avalanche  process  that  results  in  a  rapid  rise  in  the 
plasma  density  and  eventual  microwave  absorption  or  reflection.  The  continuity 
equation  for  electrons  in  the  media  is  written 

Q_ 

—  =  t/in-Uan  +  V^{Dn)  (1) 

where  n  is  the  electron  density,  Ui  is  the  ionization  rate,  is  the  attachment  rate, 
and  D  is  the  electron  diffusion  coefficient.  Usually  the  diffusion  term  is  approxi¬ 
mated  as  V^(Dn)  =  n(—D/A^),  where  A  is  the  chairacteristic  diffusion  length  and 
is  dependent  upon  the  particular  system  geometry.  Under  this  approximation,  the 
continuity  equation  yields 


n^e 


ut 


(2) 


where  i/  is  a  net  rate  for  electron  production.  A  reasonable  estimate  of  the  time 
required  to  reach  breakdown  is  obtained  if  it  is  assumed  that  breakdown  occurs 
when  the  plasma  density  reaches  the  critical  plasma  density  for  wave  propagation 

rip  =  (cm”^)  (3) 

The  pulse  duration  r  for  breakdown  to  occur  is  then  given  by 

Hp  =  ri^c"*  (4) 

Although  the  initial  free  electron  density  is  not  known  to  any  degree  of  accuracy  in 
most  experiments,  the  exponential  nature  of  the  rise  in  electron  density  with  time 
tends  to  minimize  the  effects  of  errors  in  the  number  assiuned  for  Ug.  Nevertheless, 
this  effect  can  be  significant.  For  3  cm  microwaves,  for  example,  the  difference 
between  rio  =  1  and  rig  =  1000  cm"^  represents  about  a  25%  change  in  the  break¬ 
down  pulse  duration  assuming  other  quantities  remain  fixed.  Of  course,  if  a  media 
can  be  found  in  which  rig  —  0,  breakdown  will  not  occur  for  any  pulse  length  un¬ 
less  the  wave  electric  field  is  large  enough  to  directly  strip  an  electron  from  a  neutral. 
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III.  Experiments 

High  power  microwave  pulses  used  for  the  propagation  experiments  were  gener¬ 
ated  by  a  large-orbit  gyrotron  operating  at  10  GHz,  as  reported  previously®.  This 
device  produces  high  power  microwave  radiation  (100-500  MW)  by  the  resonant  in¬ 
teraction  of  a  rotating  electron  beam  with  the  modes  of  a  muitiresonator  magnetron 
circuit.  The  experimental  configuration  used  for  the  propagation  experiments  is 
shown  in  Fig.  1.  Microwaves  from  the  large-orbit  gyrotron  aie  refiected  90®  by  a 
plane  refiector  toward  a  150  cm  diameter  parabolic  reflector.  The  parabolic  reflector 
serves  to  focus  the  microwaves  into  the  test  cell.  The  test  cell,  a  30  cm  diameter 
cylindrical  acrylic  vacuum  chcunber,  is  fitted  with  a  high  accuracy  vacuum  gauge 
and  a  calibrated  leak  valve  that  allowed  maintenance  of  ambient  pressures  from 
atmosphere  down  to  0.1  Torr.  Measurements  of  the  microwave  power  density  along 
the  dish  axis  (Fig.  2)  and  across  a  cell  diameter  (Fig.  3)  indicate  that  a  reasonable 
focus  was  achieved  away  from  the  test  cell  walls  (surface  breakdown  on  the  test  cell 
wadis  must  be  avoided  if  true  media  breakdown  limits  are  to  be  established). 

The  shield  wall  between  the  lairge-orbit  gyrotron  and  the  test  cell  was  constructed 
of  20  cm  thick  high-dcnsity  concrete  block  augmented  by  1  cm  of  lead  sheet.  Mea¬ 
surements  of  the  x-radiation  observed  at  the  test  cell  indicate  that  the  shield  wall 
reduced  x-radiation  by  more  than  two  orders  of  magnitude,  with  typicad  TLD  de¬ 
tector  dosage  recordings  of  50  mR/shot  observed  without  the  shield  wall  and  0.5 
MR/shot  after  the  shield  wadi  was  constructed.  Because  of  reflections,  it  wais  not 
possible  to  completely  eliminate  all  x-radiation  from  the  test  cell  area. 

Open  ended  X-Band  (8-12  GHz)  waveguide  sections  were  used  to  detect  the 
microwave  signals  before  and  adter  the  focal  point  in  the  test  cell,  ais  shown  in 
Fig.  4.  Microwave  signals  were  passed  through  calibrated  attenuators  to  cadibrated 
J.-loctors  oclt-clcu  for  their  measured  trainsient  response  time  (rise  time  ~  1  ns). 
Detector  signads  were  fed  directly  into  Tektronix  7104  oscilloscopes. 

In  this  work,  brcaJcdown  of  the  media  was  said  to  occur  if  either  visible  light 
from  ionization  processes  was  observed  in  time-integrated  photographs  of  the  test 
cell  or  if  significant  shortening  of  the  post-focus  microwave  pulse  was  observed  due 
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to  reflection  and/or  absorption  of  the  microwaves  by  the  breakdown  plasma.  The 
criteria  for  pxilse  shortening  to  indicate  breakdown  was  if  the  trzmsmitted  pulse  du¬ 
ration  wzis  80%  or  less  of  the  incident  pulse  duration.  Microwave  power  density  was 
varied  over  the  range  50-300  kW/cm^  and  pulse  duration  was  varied  over  the  range 
3-20  ns.  Data  for  the  experiments  axe  plotted  in  the  familiar  E/p  vs  pr  format  in 
Fig.  5. 

rV.  Conclusions 

From  the  data  shown  in  Fig.  4  and  previously  reported  data®,  it  appears  that 
a  significant  reduction  in  background  radiation  does  increase  the  wave  electric  field 
required  for  breakdown,  at  least  at  high  (>  1000)  and  low  (<  10)  values  of  pr. 
This  increase  in  the  breaJcdown  electric  field,  about  a  factor  of  ;..vo  over  the  previ¬ 
ously  reported  data,  is  small  but  nevertheless  large  enough  to  ensure  that  statisticcd 
fluctuations  aJone  cannot  account  for  the  increeise. 

Unforttmately,  it  is  nearly  impossible  to  measure  accurately  the  actual  free  elec¬ 
tron  density  in  the  test  cell  at  the  time  these  measurements  were  made,  due  pri¬ 
marily  to  the  unavailability  of  diagnostics  sensitive  enough  to  yield  quantitative 
estimates  of  such  low  electron  densities.  Even  a  theoretical  estimate  would  require 
a  very  careful  measurement  of  the  x-ray  energy  spectrum  zuid  information  on  x-ray 
induced  ionization  cross  sections  over  the  entire  spectrum.  Nevertheless,  the  data, 
and  the  theoretical  discussion  in  section  II,  suggests  that  breakdown  thresholds  as 
reported  here  and  in  the  previous  report*  are  reaisonably  reliable  as  long  ais  at  least 
a  few  free  electrons  are  present  in  the  media.  Natural  radioactivity  and  cosmic 
radiation  are  probably  sufficient  to  ensure  that  this  is  the  ceise. 
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Figure  Captions 


1.  Experimental  configuration  used  for  shielded-source  microwave  propagation 
experiments. 

2.  Schematic  of  microwave  diagnostic  setup  used  in  the  propagation  experiments. 

3.  Microwave  power  density  vs  distance  along  the  axis  of  the  parabolic  refiector. 

4.  Microwave  power  density  vs  distance  across  a  diameter  of  the  test  cell. 

5.  Microwave  propagation  data. 
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